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Photograph taken C-I-L Central Research Laboratory, McMasterville, 


Plastics Detective 


This photograph shows the High Shear Viscometer, designed and built our 


Central Research Laboratory, used investigate factors involved complex 
flow behaviour molten thermo- plastics during commercial processing. 
Information this newly developed pressure-driven capillary Viscometer 
has been made available industry C-I-L, and many industrial laboratories 
the United States, England, India and Australia have built identica! machines. 


Such C-I-L research contributes higher standards performance the 
plastics industry, new and better plastics products and—through better technical 
service—to greater success the application these new materials modern living. 
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The Jane 


Memorial Lecture Award 


Terms Reference 


The Jane Memorial Lecture Award was established commemorate 
the memory the late Dr. Robert Stephen Jane who has made outstanding con- 
tribution the chemical profession and the chemical industry Canada. 


The Award shall presented resident Canada for exceptional achieve- 
ment field chemical engineering industrial chemistry. 


Award shall presented the Annual Conference The Chemical 
Institute Canada. The recipient the Award shall required present lecture 
the Annual Meeting the Canadian Chemical Engineering Conference The 
Chemical Institute Canada the year which the Award made. The subject 
the lecture shall the outstanding work the lecturer the field chemical 
engineering industrial chemistry. The lecture will published The Canadian 
Journal Chemical Engineering. 

candidate chosen for this honor shall known “The Jane Mem- 
orial Lecturer”. will receive cash award $300, and travelling expenses. 

intended that the Award shall granted annually. However, order 
maintain the standard the Award, the Selection Committee may omit the Award 
any year which suitable recipient not available. 

The recipient may become eligible for the Award subsequent year, pro- 
vided that has again made new significant contributions chemical engineering 
industrial chemistry Canada. 

Nominations for the Award shall made writing the Secretary The 
Chemical Institute Canada, Rideau St., Ottawa Ont., not less than five (5) 
professional members the Institute, January each year. All accepted nom- 
inations remain good standing for period three years. 

The Director Scientific Affairs shall refer all nominations the Selection 
Committee which shall consist of: 

The Chairman the Chemical Engineering Division 

The Editor The Canadian Journal Chemical Engineering 

representative the Trustees The Jane Memorial Lecture, 
selected the Trustees 

representative from the Chemical Engineering Division appointed 
the Director Scientific Affairs; the latter non-voting member 
this committee. 


The Award Committee shall determine, prior March the candidate who 
receive the Award. 


Dr. Rowzee, president and managing director, Polymer Corp. 
Limited, Sarnia, Ont. the first recipient this Award. 
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CHEMICAL ENGINEERING AWARDS 


Commencing last fall, the Chemical Engineering Division the established 
two awards the first for the best paper published The Canadian Journal Chem- 
ical Engineering during the year immediately preceding the Fall Conference, and the 
second for the best oral presentation technical paper the Annual Con- 
ference (Chemical Engineering Session) and the Fall Conference the Chemical 
Engineering Division. The awards will take the form engraved certificate and 
silver tankard with the proper inscription. 


The first award for best paper was won Torobin and Gauvin 
McGill University and the Pulp and Paper Research Institute, Montreal, for their 


series papers “Fundamental Aspects Solids-Gas Flow”. 


The first oral award winner was Smith, Northwestern University, Evanston, 
Ill. for his presentation “Heat Transfer Surface Reacting Fluid Turbulent 
Flow”. His co-author was Ronald Rothenberg, and the paper was published the 
December 1960 issue The Canadian Journal Chemical Engineering. 


CHEMICAL ENGINEERING DIVISION EXECUTIVE FOR 


Chairman: Atomic Energy Canada Limited, Chalk River, Ont. 
Vice-Chairman: McMaster University, Hamilton, Ont. 
Immediate Past Chairman: McGill University, Montreal, Que. 


Executive Members: University Alberta, Edmonton, Alta. 


Scorr, University British Columbia, Vancouver, 
Pont Canada Limited, Sarnia, Ont. 


Kemper, Dominion Tar Chemical Co. Limited, Montreal, 
Que. 


Atomic Energy Canada Limited, Chalk River, Ont. 
Craic, Ontario Paper Co. Limited, Thorold, Ont. 
Jounson, University Toronto, Toronto, Ont. 
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THE JANE MEMORIAL LECTURE 


RALPH ROWZEE, 


President and Managing Director, 
Polymer Corp. Limited, 
Sarnia, Ont. 


Rubber, Research and Human 


Resources 


honored the 1960 Jane Memorial lecturer; 
particularly since this the inaugural lecture this fitting 
memorial. the first recipient the award, aware 
heavy responsibility; namely, that presenting address 
which Dr. Jane might take pride. sincerely hope that may 
discharge this responsibility creditably, and that may conform 
the high standards which Dr. Jane always set for himself. 

knew Dr. Jane for many years, and held him high regard 
friend and fellow chemist. appreciated his quiet, unassuming 
manner, and respected his views and judgment. particular, 
should like emphasize his contribution our profession, 
our industry, and Canada. His contribution the 
testimony the faith had the chemist, the chemical engi- 
neer, and the chemical industry. 


Dr. Jane had sincere interest people and their progress. 
The status enjoyed the company which served 
reflection his ability and leadership. His thinking was always 
broad, imaginative and dynamic. His efforts extended beyond 
his company and his served his community and 
his country. should like, therefore, feel that would 
accept this address contribution and challenge the 
industry worked for and the people worked with. 

The terms reference for this lecture stipulate that the 
subject shall related the work the lecturer the field 
chemical engineering industrial chemistry. There little 
doubt therefore, that make some reference synthetic 
rubber but, may, should like interpret this memorial 
believe Dr. Jane would wish it, opportunity use the 
history the past only challenge the future. 


This Lecture was presented the Canadian Chemical 
ing Conference, Quebec, Que., Nov. 7-9, 1960. 


vii 


The title the lecture Research and Human 
should like start out, strangely enough, with 
natural rubber. not technical paper, but rather 
statement experience the rubber industry and certain 
conclusions which have drawn regarding technical education 
and research. 


Brief History Rubber 


The story rubber interesting one, repeatedly empha- 
sizing the action the economic laws supply and demand, 
and the fact that where unnatural controls are imposed, way 
will found circumvent them. 


The tale has often been told how Sir Henry Wickham, 
1876, smuggled seedlings the rubber tree out Brazil 
Kew Gardens, from where they were transplanted Ceylon 
and the Straits Settlement. Here they formed the basis for the 
far eastern plantations. 


The real stimulus which led the large scale development 
these plantations, however, was the market pool operated 
Brazilian growers and speculators which ran prices natural 
rubber over $1.00 pound 1906 and $3.00 pound 
1910. While the speculators reaped the harvest, they sounded 
the doom Brazilian rubber. The industry moved the far 
east and British interests dominated the field. 


the arbitrary price controls the British, through 
the Stevenson Act, restricted output and affected supply. The 
result this control was raise the price, encouraging heavy 
plantings other interests (primarily the Dutch). ‘It also led 
some serious thought about synthetic rubber. The immediate 
result was that independent planters became active, and more 
and more rubber came the market. The British domination 
the field declined, and the price dropped supply increased, 


With the depression the huge surpluses existed 
and rubber prices dropped low pound. agree- 
ment reached 1934 imposed restrictions international 
scale, and prices gradually recovered the 20-25¢ level. 
the last quarter 1941, just before Pearl Harbor, all rubber 
possible was being shipped out the East. This significant 
item since provided notice the urgent need for rubber that 
stimulated the wartime synthetic rubber program North 
America. 


Synthetic Rubber 


Synthetic polymers with certain rubber-like properties had 
been known for almost century, starting with the work 
Tilden the 1860’s. During World War natural rubber was 
cut off from Germany the British Blockade and the German 
scientists succeeded developing rubber-like polymers based 
butadiene which they used for truck tires and gun 
carriages. The material was not good rubber and, with the 
armistice, production lapsed until the Nazi regime came 
power with ambitions for world domination. Realizing that 
rubber would again cut off the beginning war, the full 
support the government was given the scientists evolve 
rubber-like polymer which would serve substitute for natural 
rubber. 

Among the many chemical advances made Farben 
chemists during this period, were the discoveries emulsion 
polymerization, copolymerization, and series polymers based 
butadiene. sodium catalyzed butadiene rubber was devel- 
oped. Factories were designed for copolymers butadiene and 
styrene (Buna and just before the war, were brought into 
production. 

The Standard Oil Company acquired the rights many 
the Farben patents through arrangement with Farben, 
which provided for exchange information areas mutual 
interest; thus, knowledge buna polymers came America. 

America, with plenty rubber available, there was little 
incentive for the development synthetic rubber during the 
early 30’s. The matter received attention, however, and research 
quietly went forward. 1922 Thiokol was accidentally dis- 
covered American chemist, Patrick. Thiokol was 
oil resistant, but had other less desirable characteristics, and 
although used for few special applications, never became 
important product. the late DuPont acquired the 
rights promising material discovered Father Nieuland 
Notre Dame University. This material, later called neoprene, 
excelled certain characteristics where natural rubber was 
weak, but was some years before neoprene was produced 
commercially. 


was this time, 1931, that graduated from 
chemical engineer, and went work the Research Division 
the Goodyear Tire and Rubber Company, Akron. Good- 
year had done some research work from time time synthetic 
rubber, but had always been interrupted more urgent things. 
1933, Litchfield, one the outstanding leaders 
the rubber industry, and then President Goodyear, called 
for full time major assignment synthetic rubber. Asked 
about the possibility conflict with German patents, Litchfield 
reported have said: the cross that 
bridge when come it. Our job now find out every- 
thing can about synthetic and from 


From Test Tube Test Tire 


1935, was assigned the job devising means, and then 
producing, lots approximately ten pounds each, variety 
the most promising synthetic rubbers based butadiene. 
Test tube quantities previoulsy produced Goodyear chemists 
had been insufficient for compounding evaluation. Within 
short time, the assignment was expanded the production 
several hundred pounds particular rubber that test tires 
could produced. Today, produce synthetic rubber 
such quantities that difficult visualize the time and work 


involved producing few pounds rubber 1936, let alone 
several hundred pounds. 


Some the problems faced those days were: long reactions 
often requiring weeks complete; reactions were stopped short 
completion because accident was discovered that 
70% conversion monomers tough rubber-like material was 
produced whereas 90% conversion, the material had few 
rubber-like properties. Reasons for this were 
understood the organic chemist but were later explained 
the physical chemist; there were facilities recover 
unreacted monomers; equipment had improvised for 
coagulating, washing and drying the product. 

Other problems presented major assignments themselves. 
the principal monomers, only styrene could purchased. 
Butadiene, acrylonitrile and methylacrylonitrile had pro- 
duced from intermediate materials the laboratory, using very 
crude and laborious techniques. 


Many years after the event, engineer the 
Company Elyria, Ohio, recalled first visit his plant 
discuss reactor for synthetic rubber. seems that repre- 
sented one two Akron rubber companies who approached 
him within few days one another. Each presented sketch 
five gallon glass-lined pressure vessel which Pfaudler was 
manufacture quickly possible. said that was 
amazed find the sketches similar. Those five gallon 
reactors produced Pfaudler 1936 were the first what 
later became standard item the laboratory every copolymer 
plant. From these developed the large, glass-lined reactors which 
are found copolymer plants today. 

This basic development period was followed period 
scaling through pilot plant, and the construction two 
plants; one that Goodyear built for the United States Govern- 
ment which was six months away from completion the time 
Pearl Harbor, and one that Goodyear operated privately. 
retrospect, stands out clearly that most the things that had 
done this early development stage, made use 
chemical engineering education. 


Date with History 


Shortly after Pearl Harbor, the synthetic rubber program 
was adopted and urged with high priority. spite the 
plentiful supply natural rubber available few years earlier, 
North America suddenly became continent the 
matter supply this important commodity. These words 
“have were used the Rubber Study Committee chaired 
Bernard Baruch 1942 stress the urgency the situation 
with which were faced. The Committee, its report, 
emphasized that large scale production synthetic rubber was 
the end 1943. Victory defeat was the balance. 


The Canadian Government decided early 1942, par- 
ticipate this program one its major contributions the 
war effort. undertook build plant Sarnia produce 
both GR-S and rubbers, well butadiene, isobutylene, 
and styrene, the major chemicals required the processes. 
was assigned the copolymer part the project, loan 
from Goodyear. This was the start association with 
Polymer, and with Canadians, and the beginning becom- 
ing Canadian. 

The synthetic rubber program has been called one the 
major chemical engineering achievements all time. Certainly, 
was not without its problems that now provide interesting, 
pleasant reminiscence. major significance was the coordina- 
tion pooled processes, and ideas many companies, into 
so-called standard plant. The development the standard plant 
excellent example good chemical engineering design. 
more remarkable because the limited operating experience 
that had been gained before Pearl Harbor. 

How good was the have been many pieces 
equipment that, for one reason another, operators and 
management would have liked replace, but spite real 
effort find better equipment, relatively few changes have 
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made. Most basic equipment ordered for new plants to- 
very similar the equipment installed existing plants 
almost twenty years ago. 

This does not imply, course, that changes have not taken 
place. the contrary, many improvements have been made, 
but these have been primarily process changes, and changes 
auxiliary equipment. The conversion the copolymer reaction 
from batch continuous process was highly important 
since increased the capacity the plant substantially. New 
controls and instrumentation have improved uniformity and 
other quality characteristics; new types rubber have been 
introduced supply the specialized needs industry. These 
are but few the changes, and all them can seen the 
work the chemical engineer. 


Polymer has pioneered many new types rubber, and has 
made many noteworthy industrial chemical and chemical engi- 
neering achievements which give the Company important tech- 
nical advantages, some which have not yet been disclosed. 
The application Dow, type dehydrogenation catalyst 
example the constructive work done this field the 
technical staff the Company. 

Further indication the application chemical engineering 
work has been the increase capacity achieved bottleneck 
removal programs and process improvements. 


The Construction Period 


should like refer briefly the construction period. 
The first rubber was produced September 29, 1943, well 
ahead the deadline established the Baruch Committee 
its report, which referred earlier. That were able 
tribute all who participated; engineers, contractors, 
suppliers, administrators and workers all kinds. true, 
there was great incentive, but the feat was accomplished 
the face great odds and we, Canada, can justly proud 
our participation. 

Many people considered the synthetic rubber plants 
expendable items wartime emergency, and that recovery 
the natural rubber areas would make these plants orthless 
least idle and could have been. The challenge, 
however, was there make industry for peace-time. 
meet this challenge, was essential find market for the 
capacity the Polymer plant; capacity that could supply 
much more than the Canadian requirements for rubber. 


Development Markets 


The places turn, course, were the major rubber con- 
suming areas the world. proposal develop these export 
markets was accepted, though most advice said that Europe 
would never buy synthetic rubber. Starting Europe, top 
Polymer representatives dealt with trade commissioners, govern- 
ments and rubber consumers. Their efforts were rew arded and, 
eventually, orders started flow. Today, the western world 
consumes almost 4,000,000 tons rubber annually, compared 
with 1,000,000 tons years ago and almost half this present 
synthetic. Polymer exports 65- 70% its 
output, and its products flow more than countries the 
This annual export product valued excess 
$50 million, and important factor Canada’s economy. 

This has required concentrated effort marketing and 
technical service. has required consistency approach 
supply and price. has required diversification production 
and strong research effort. Starting 1946, competent agents 
were appointed major countries, and this group has made 
significant contribution the respect which the name 
POLYMER held. 

this highly competitive industry, Polymer operates 
manner similar that any private company. pays taxes 
and dividends, and absorbs all the normal costs business. 
contributes its full share the life and welfare the commu- 
nity which operates. the last fifteen years, using the 


same basic facilities, production has been increased from 50,000 
150,000 tons per year. 


Perhaps not too much say that this amazing 
record. 1941, Canada was most unlikely location for 
synthetic rubber plant, and one would have been considered 
out his mind had predicted that within fifteen years 
Canada would significant factor the world rubber supply. 


believe the story Polymer points out the possibilities 
enterprise against heavy odds. highlights the contribution 
the chemist and the chemical engineer meet the challenge 
“have situation with great benefit the 
country. can said that the circumstances were unusual, 


but ideal conditions are most unlikely exist for any venture. 


Canada growing nation. Estimates population 
over million 1980 have been given comprehensive 
studies Canada’s future. the present time, Canada has 
wealth natural raw materials, but because our relatively 
small population, considered too small market for the 
mass production conversion out natural resources consumer 
goods. This is, sure, sound economics consider the 
limits the Canadian market but need look 
such? Are there not products that can produce for 
much larger market, and should not look wider fields 
sell products which could be, but have not yet been, produced 


Canada? 

Our prospects for prosperous Canada, with population 
million, are dependent the opportunities can provide 
for employment, not only our present industries, but also 
new industries and services which should begin take shape 
the next ten years. This growth will bring Canadian 
market better able sustain many the consumer goods 
industries, but there bridge cross before this achieved. 
cross will require the utmost initiative and enterprise 
all Canadians; management and labour financier and 
investor manufacturer and consumer educator and in- 
dustry government and the people. The challenge great 
are prosper generally buoyant world economy. 


Important Factors 

not have any panacea assure success this growth. 
know that will mean hard work, understanding and coopera- 
tion among all segments our people. should like consider 
two factors that shall need the next decade, which 
believe will make special contribution the development 


Canada. These important factors are: 
One our trained people, and 
Two industrial research. 


think that will evident you that these two fields are 
particular interest because industrial association 
and background. 


First, let discuss our trained people! 


Since our economy based the second highest average 
labor cost the world, technology great importance. 
order compete with goods manufactured lower cost labor, 
must concentrate discovering, ahead other countries, 
know-how many appropriate fields possible. must 
also strive for the highest productivity level can achieve 
automation, work simplification, and plain hard 

Canada already owes much the chemist and the chemical 
engineer. Even more will depend these people the future. 
You are present this meeting because your interest 
technical development. You have all graduated from univ ersity; 
most you from Canadian universities, and many you 
recently. should like you consider with some aspects 
this educational system. 

have, Canada, excellent educational system all 
levels. There will always arguments about speeding the 
progress the above erage student through elementary and 
secondary schools. There the problem facilities meet 


; 


the increasing enrollment students. Industry has recognized 
the need the universities for expansion turn out more 
trained people and, along with other groups, has contributed 
generously the building program our universities. 


What About Our Courses? 


What has been happening the curricula our engineering 
courses during the last several years? believe that the technical 
training has many cases been diluted what will call 
secondary factors. There has been tendency for employers 
demand, and for the universities adopt, courses which 
provide such generalized training that the individual graduating 
from these courses lacks depth technical training fill the 
engineering chemical requirements industry. Much 
this change has been introduced with the objective graduating 
managers. There are many people who believe this possible. 
not disagree with the desirability producing managers 
from educational system but, believe that are years 
away from perfecting the techniques accomplish this. 
prime importance such development the selection the 
right candidate. the meantime, these generalized courses 
tend develop impatience the graduate become execu- 
tive over night. You and know that graduate cannot become 
effective supervisor, let alone executive, until has met 
the tests that only experience can give him. 


are all blame for this development. The universities 
have been under heavy pressure from industry 
generally, provide such training. suggest that good 
scientific training, specializing the field the individual’s 
particular ability, prime importance the engineer and 
the chemist. For those especially talented leaders managers, 
further training the graduate level university, within 
industry, will better equip them meet the challenge than will 
the generalized type undergraduate course. Some authorities 
are now pointing out that additional training graduate level 
would most effective after some years experience 
industry. 

These remarks have been directed toward under- 
graduate courses chemistry and chemical engineering. 
must not, however, lose sight the humanities. well 
aware the time factor, but believe that the graduating 
student who has had introduction the humanities should 
better adjusted meet the challenge his profession. 
would make clear that believe emphasis must first placed 
the technical aspects the course, but that the student 
should not permitted lose sight the importance the 
humanities, particularly with the continually broadening sphere 
activity and employment available the chemical engineer 
and science graduate. 


Compensation 


Another important development which taking place relates 
graduate employee performing top flight technical work 
accord with salaries paid for supervisional and managerial 

work. This, believe, properly takes away some the 
anxiety the young graduate who seeks the utmost financial 
reward through early advancement managerial posts, and 
provides better perspective the true importance the 
scientist. There delicate balance achieved this 
matter compensation, and that the scale now 
tipping achieve this balance more appropriately than has been 
the case past years. 


When speak productivity, are too apt think 
only the output our wage earners. Excellent articles have 
been written present the argument that productivity 
measure the total output achieved from total 
important that should not lose sight the fact that manage- 
ment, research, purchasing, traffic, accounting, and 
marketing along with production, fact, all costs business 
are part input contributing the productivity the 


business. therefore essential that utilize all 
effectively achieve the utmost productivity. 


mention this because there has been tendency use 
chemists and chemical engineers the capacity technicians, 
which hardly keeping with the objectives using people 
most effectively. believe mistake have “chemist 
chemical engineer doing routine work that could done 
regularly less qualified individual. Even the chemist, 
chemical engineer, prepared accept such assignments, 
misuse the training developed our educational 
therefore desirable stress the proper utilization our 
trained people industry and business. 


Today the chemist chemical engineer being used 
Management and Marketing, Economic Evaluation and 
Technical Services, and Research and certain special functions; 
being used wherever his particular talents indicate that 
there activity through which can make significant 
contribution his organization. This represents great change 
from the situation which prevailed when entered industry and 
that which found coming Canada almost years ago, 
when was told number people that the chemical engineer 
was neither fish nor fowl, that was neither chemist nor 
engineer, and therefore little value. 


The Use Trained People 


There one other aspect the use our trained people 
which would refer. have, over the years, tended 
confine the chemist the laboratory and the chemical engineer 
the plant. This probably results from looking the chemical 
engineer the applied chemist and the plant problems applied 
problems. Much might gained assigning the chemist 
the plant and the chemical engineer the laboratory, not 
necessarily try interchange the two individuals, but rather 
bring fresh perspective the problems the plant and the 
problems the laboratory. Improvements processes and 
process design might result from the fresh approach the 
chemist and improvements laboratory equipment might 
developed the engineer the laboratory. Moreover, 
having both chemist and engineer working together cooperatively 
and objectively provides combination talent that might 
recognize new concepts and important departures from our 
standard processes that will contribute the industrial growth 
seek. 


The Case for Research 


Research the second factor that should like have you 
consider. You are all well aware the importance 
research technological development. Considerable research 
has been carried Canada the universities, hospitals, 
research institutes, industry, and especially the National 
Research Council. research, with which 
ticularly concerned, has stimulated interest specific fields 
research and this has led, turn, the stimulation ‘basic 
research related fields within the universities. example, 
believe that research carried Polymer Corporation has 
stimulated interest polymerization generally, which turn 
has had bearing the the High Polymer Forum 
and the Polymer Research Institute. This healthy kind 
development. 


Industrial research any company, course, justified 
the basis that will provide return that company. 
There is, unquestionably, risk involved. The stakes are great, 
and the payout can most attractive. Polymer has benefited 
from the research carried its staff, both the technical 
and process fields, and the respect which accorded its 
research and technical people. Polymer has 
from the research being carried the N.R.C. the univer- 
sities, and from the work the High Polymer Forum. 

should like stress, however, that with Canada’s economy 
tied closely that technological development, research 
ascending scale important, not only the interests 
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extracting raw materials, but also converting our raw materials 
consumer goods. 

The case for industrial research has been advanced before. 
should like review some the arguments supporting this case: 


Research required develop know-how and world-wide 
patent positions. 


New discoveries that lead knowledge and experience 
any field advance other countries and, that can upheld 
strong patent position, have far reaching effects. They 
stimulate further research the same and related areas. 
addition the financial reward, the resulting prestige adds 
strength negotiations with others and leads exchanges 
with others doing important work the same basic fields. 


Research required develop new equipment, new ma- 
chines and new techniques. 


have mentioned the importance productivity the field 
competitive marketing. must not satisfied continue 
doing things the same way, with the same tools, but must 
continually seek for better ways and means. Our ingenuity 
developing machines and techniques may contribute greatly 
the productivity factor. This primarily job for industrial 
research. 


Research required develop secondary and consumer 
goods manufacturing Canada. 


This development closely related the size the market. 
long accept the idea that the market can develop 
Canada and other areas not large enough for operation 
Canada, can expect dependent other countries 
and, therefore, provide employment opportunities those other 
countries rather than Canada. Industrial Research one 
the many factors that can contribute the development 
secondary and consumer goods manufacturing. 


Research required encourage the interest young 
people scientific careers and, through enhanced stature, 
encourage our talented, technically trained people remain 
Canada. 


have emphasized the excellent educational system Canada, 
and that science and engineering graduates from this system 
play important role Canada’s future. industrial research 
assumes broader and more important role, the stature those 
engaged will greatly enhanced. This important 
factor attracting young people ability, and the right 
temperament, into work. 

These arguments are interest the Canadian owned 
firm. Are they equal interest the company whose Canadian 
plant subsidiary foreign controlled corporation which 
conducts centrally operated research facilities its home terri- 
tory? This is, course, difficult situation. Every owner 
must weigh the pros and cons, and determine the best course 
action for his own interest. suggest that should 
continually aware the contribution which research Canada 
can make this country with material benefits himself when 
planned and conducted the basis long term approach. 

Several important Canadian companies whose control lies 
outside Canada have taken this approach, believe their 
advantage. 

have mentioned that Dr. Jane was not restricted his 
thinking. This was particularly true his views research. 
recall discussing with him ways and means stimulating re- 
search Canada. was prepared explore every avenue 
which might have chance success. advanced the idea 
that several companies might pool their research activities for 
the benefit all, and Canada. not intention discuss 
thought this kind. However, would bring your attention 
the need for imaginative approach Canadian needs and the 
means fulfilling ‘them; approach that may generate the 
growth Canada industries now have, and new ones 
that may appear uncertain success Polymer did the 
early 

believe Dr. Jane, all his many activities, was thinking 
and working for better Canada. worked for the develop- 
ment and growth industry and society. hope that 
remarks may accepted the same spirit, and that they would 
acceptable him. 
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new type mixer-settler described wherein 
the heavy liquid phase spouted upward through 
depth light liquid phase and allowed settle 
below the light phase. The method does not require 
mechanical agitation, allows both mixing and settling 
continuously the same chamber, and 
simple design. 

Optimum conditions have been determined for 
extraction acid from toluene 
means water spout. Independent variables 
studied were water rate, toluene rate, spout pipe 
diameter and wall effect. Dependent variables meas- 
ured were spout height, spout spread, initial entrain- 
ment points and stage efficiency. Murphree spouted- 
phase efficiencies 96% were obtained. 


single-stage mixer-settlers, when they are con- 
tinuous flow devices, almost invariably require that the 
mixing and the settling separate units 
This because the mixing brought about 
either mechanical agitation parallel flow the 
two liquids being 

The present device was originally conceived two-pass 
cocurrent-countercurrent liquid-liquid contactor. was sug- 
gested and named after the fluid-solid contacting operation 
developed Mathur and though the resemblance 
between fluo-solid and the present liquid-liquid 
operation superficial. The intention was inject the denser 
liquid high velocity vertically upward through cocurrently 
flowing column the lighter liquid, the interface which with 
the heavier liquid was located immediately below the 
heavy liquid inlet. was hoped that the spout dense liquid 
would rise through very large height light liquid before 
breaking into droplets, which would then descend counter- 
currently through the annular region the upward moving 
light phase. this way two long passes the heavy phase 
were perform the equivalent several mass transfer equilib- 
rium contacts. 

fact, however, the liquid-liquid spout heights high 
velocity were considerably smaller than, for example, those 
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obtained spouting liquid into air (Figure spouting 
velocities were limited entrainment considerations, and 
the case confined column spouting (Figure Ic), turbulently 
mixed dispersion zone was set the region the spout 
and below it, rather than orderly cocurrent-countercurrent con- 
tacting the two phases. result the equivalent single 
equilibrium contact was never exceeded. For this reason the 
method was studied with the altered viewpoint optimizing 
its performance single-stage single-chamber continuous flow 
mixer-settler. 


Spouting Behavior 


That liquid spouting feasible was first demonstrated 
injecting water into both gasoline and toluene via centrally 
located vertical mm. tube 2.3-in. glass column with 
contracting conical bottom, geometrically similar fluo-solid 
spouting column. was quickly realized, however, that the 
conical bottom fulfilled useful function the liquid contacting 
operation and was actually detrimental the settling and separa- 
tion the denser phase. The cone was therefore eliminated 
the subsequent column experiments, and instead, tapered 
expanding section was used for decelerating the outlet water. 


Figure schematic diagram the apparatus used for 
studying the characteristics water spout toluene 
4-in. Pyrex column under both continuous-flow and semi- 
batch (zero toluene flow) conditions, and inverted 13-gal. 
glass carboy under semi-batch conditions alone. Except for the 
copper breather tube the carboy, which caused slight dis- 
coloration the toluene with time, especially the presence 
dissolved benzoic acid, and some ace Saran sampling leads, 
all lines were type 316 type stainless steel. standard 
tapered 6-in. 4-in. Pyrex reducer the water outlet was 
used minimize bottom entrainment toluene water 
the column, while 6-in. toluene outlet section was employed 
diminish any overhead losses water. The narrow 3-in. 
diameter neck the bottom the 16-in. diameter carboy, 
the other hand, served enhance bottom entrainment toluene, 
thereby limiting the carboy water capacity 1.2 
gal./min. 

Interface control, initially effected barometric loop 
from the water outlet, was later accomplished manipulating 
the outlet water valve. Where well-defined toluene-water 
interface was formed, the carboy runs and the beginning 
column run, this interface was kept slightly below the top 
the spout pipe prevent its disturbance the high 
velocity water leaving the pipe. For the continuous ‘flow runs, 
toluene was introduced the initial interface level via two 
1/4-in. schedule pipes adjacent the spout 
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Figure representation water spouting 
into (a) air, (b) unconfined toluene, (c) confined toluene. 


facial scum was eliminated between runs, was any discolored 
toluene the carboy. 

Figure shows spout height function water rate, 
zero toluene flow the 4-in. column, for o.d. and 
1/8-in., 1/4-in., 3/8-in. and 1/2-in. schedule spout pipes. 
For pipe, with the possible exception the largest, spout 
height rises approximately linear fashion first and then 
tapers off with water rate. given volumetric water rate, 
spout height increases with decreasing inside diameter the 
spout pipe, would expected from the increase 
entering velocity the spout. low water rates spouting 
poor there fountain-formation the peak the 
ascending water column, which breaks into large globules 
rather than dispersing into small drops. This region 
followed one good spouting, which well-formed 
fountain the summit the spout breaks into well-dispersed 
droplets. Finally limitation the region stable spouting 
reached when some the water droplets formed are 
small entrained upward the slightest upward flow 
toluene. should noted that the water flow range 
which good stable spouting occurs, indicated the region 
between the dotted lines, increases with spout pipe size first 
and then decreases sharply. 


The effect upward toluene flow rate spout height for 
1/8-in. schedule spout pipe shown Figure for 
two constant water rates. The small influence toluene rate 
can explained the fact that the average linear velocity 
the toluene only 1/220th that the entering water spout 
velocity, for equal volumetric rates toluene and water. This 
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Figure 3—Water spout heights static toluene within 4-in. 
column. 


Figure 2—Schematic flow diagram apparatus. Dotted 
line disconnected and toluene-phase supply shut off during 
semi-batch run. Tapwater supply line; Toluene-phase 
feed tank with pump; Toluene-phase receiver; 
phase drain lines; Water-phase sampling lines; 
Toluene-phase sampling line; Breather pipe; Rota- 
meters; Dispersion-water interface; 10. Air-toluene inter- 
face; 11. Toluene-water interface. 


explanation could also account for the larger effect toluene 
rate the lower water rate, where the linear velocity toluene 
relative the water spout greater factor two than 
the higher water rate, for given volumetric rate toluene. 


The glass carboy was used observe under semi-batch 
conditions the characteristics unconfined spout, diagram- 
matically illustrated Figure 1b. Here the water spout breaks 
into shower discrete droplets which flare out they descend 
into the relatively quiescent toluene medium. The photographs 
Figure give some impression how globule-formation 
superceded good spouting water rates exceeding 0.5 
gal./min., and how the resulting fountain continues spread 
out radially with increasing water rate, the same time its 
vertical growth arrested. These trends are also shown 
Figure where radial spread and vertical height the uncon- 
fined spout are plotted function water rate. 

Figure also indicates that spout heights the 4-in. column 
were somewhat lower than corresponding spout heights the 
carboy, for the higher water rates. This effect believed 
caused the greater strength and proximity the spout 
the downflow currents, created the water drops falling 
within the confined annular space the column, than the corre- 
sponding currents the carboy, where the water drops are 
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Figure 4—Effect toluene rate water spout heights 
4-in. column, using schedule-40 spout pipe. 


The Canadian Journal Chemical Engineering, February, 1961 


Figure 


free 
dispersic 
and 
also accc 
water 
carboy 


The 
energy 
spout 
peak 


For wat 
spout 
low side 
spout 
(Uay_/U 
sin 
other 
are far 
believed 
consider 
absent 
with the 
the poss 
and the 


Mass 


For 
system 


The 


0 10 — 
} 
g 
1.0 
a“ 


1.0 U.S. GAL./MIN. 1.2 U.S. GAL./MIN. 


Figure 5—Photographs unconfined spout. Rates are for 
water spouting into static toluene via schedule-40 
spout pipe. 


free spread out. These downflow currents toluene 
dispersion exert retarding shear the ascending water spout, 
and sweeping with them the smallest water droplets, could 
also account for the observed fact that overhead entrainment 
water droplets occurred slightly lower water rates the 
carboy than the column. 

The dotted lines Figure represent heights spout 
predicted the following frictionless form the mechanical 
energy balance, based the space-average velocity leaving the 
spout pipe and derived assuming zero water velocity the 
peak the spout: 


For water spouting into air, Equation (1), using values the 
kinetic-energy factor computed and reproduced 
McCabe and gives reasonably good prediction 
spout heights, though the predicted values are slightly the 
low side. Had been based the axial velocity leaving the 
spout pipe, equivalent taking Equation (1) equal 
then the predicted values would have been high 
similar degree. For water spouting into toluene, the 
other hand, the theoretical spout heights given Equation (1) 
are far above the observed values. This large discrepancy 
believed caused the energy consumed creating the 
considerable new droplet surface the head the spout, 
especially the higher velocities, effect which almost 
absent during the upflow water into air. 

the large inlet velocity the base the spout, coupled 
with the large surface created the summit, which makes for 
the possibility good mass transfer between the spouted liquid 
and the surrounding liquid. 


Mass Transfer 


For studying mass transfer liquid-liquid spouting, the 
system water-benzoic acid-toluene was chosen. 


This system, 
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Figure 6—Water spouting characteristics static media 
for schedule-40 spout pipe. 


the equilibrium data for which are readily available and show 
little effect benzoic acid concentration the immiscibility 
water and toluene, has been used many other investi- 
standard for evaluating liquid-liquid extractors. 

The direction extraction was from toluene the water 
spout, which entered free benzoic acid. Four continuous 
flow runs were performed the column, with the toluene phase 
entering point slightly below the top the spout pipe. Ten 
semi- -batch runs were made, two which were done the 
carboy. Periodic sampling inlet toluene, and outlet toluene 
and water streams, presented problem the continuous runs, 
but the semi-batch runs, special lines had provided for 
simultaneous sampling both changing phases, shown 
Figure toluene phase the runs was sampled 
before the start run and quickly flushing out the carboy 
the end run. 

Water streams were analyzed for benzoic acid alkaline 
titration phenolphthalein end point. Benzoic acid analysis 
toluene streams was effected either titration with excess 
aqueous sodium hydroxide followed back-titration with 
hydrochloric acid, after Mayfield and direct 
titration with alcoholic sodium hydroxide, after Appel and 

The index used for evaluating and comparing mass transfer 
efficiency was the Murphree extract efficiency illus- 
trated Figure the continuous runs, the Murphree 
raffinate efficiency and the more fundamental fractional 
stage efficiency could also evaluated, but was far 
more sensitive experimental error and fluctuations from 
steady state than the other two efficiencies, which were practi- 
cally equal each other all cases. For the semi-batch runs, 
only could evaluated, using the simplified formula for 
the present case solute-free inlet water stream: 


Since, the extent that deviated from all, was 
always slightly smaller than for the given operating system, 
the reported results are not slanted favor the equipment. 

The equilibrium line 59.9°F. Figure that Row, 
Koffolt and who also provided equilibrium data 
four other temperatures. interpolate these data for other 
temperatures factor was defined 


where both values apply the same value 
analysis all the equilibrium data provided Row al. 
indicated that was independent and therefore 
function temperature only. This function plotted Figure 
necessary only read off from the equilib- 


from Figure 
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Figure Equilibrium line, operating line, and stage 
efficiencies for continuous run. 


compilation all the mass transfer results given 
Table Although runs and were performed the 
water rate which overhead entrainment water droplets 
toluene starts, water droplets were actually carried out 
the column continuously flowing toluene, due the quiescent 
expanded section the top. 

For the continuous runs, the first column outlet toluene 
concentrations Table are average measured values, while 
the second column are values calculated material balance 
from knowledge the steadier concentrations the other 
three streams and the relative flow rates the two phases. 
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Figure 8—Temperature correction factor for equilibrium 
data. 


The values Murphree extract efficiency based 
two figures are satisfactory agreement and were therefore 
averaged. 

For the two carboy runs, the first values and 
listed are based samplings the beginning 5-minute 
run and are unambiguously reliable. The second bracketed 
values, however, were based samples the mixed toluene 
phase the end run, and were calculated the assumption 
that all the toluene within the carboy participated the mass 
transfer benzoic acid the water. both cases the bracketed 
values are lower than the initial values. The final data were 


TABLE 
SUMMARY TRANSFER RUNS 
column 1.5 0.5 0.653 12.21 10.92 0.717 91.1 92.7 
column 2.5 0.5 0.191 12.20 11.68 0.721 26.5 26.8 
11.24 0.707 27.0 
9.78 0.656 55.9 
column 0.269 1.0 0.56 0.345 11.19 10.04 0.666 51.8 51.1 
0.279 (10.20) (0.675) (41.3) 
9.80 0.661 42.2 
carboy 0.5 0.314 9.10 9.10 0.704 44.6 
0.250 (8.96) (0.699) (35.8) 
5.93 0.557 44.9 
column 0.180 0.6 0.387 8.33 8.33 0.572 67.6 
0.314 6.87 0.518 60.6 
column 1.0 0.378 8.82 8.82 0.625 60.5 
0.357 7.93 0.540 66.1 
0.327 7.08 0.508 64.4 
column 1.5 0.488 8.50 6.01 0.507 96.2 
column 0.269 0.8 0.198 8.82 8.82 0.618 32.0 
0.165 8.53 0.606 27.2 
0.137 8.02 0.593 23.1 
0.414 22.8 43.0 
0.405 21.4 0.933 43.4 
column 1.6 0.168 10.89 10.89 0.709 23.7 
0.160 10.64 0.700 22.9 
0.110 9.46 0.619 17.8 
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Figure 9—Effect toluene rate stage efficiency 4-in. 
column, using schedule-40 spout pipe. 


therefore recalculated material balance the assumption 
that only the toluene occupying cylindrical volume length 
and diameter equal that the spout height and spread 
respectively (see Figure 1b) participated the mass transfer 
significantly during the brief interval run. The excellent 
agreement the resulting values listed immediately 
below the bracketed values Table with their respective 
initial values, appears vindicate this assumption. 


The same problem was avoided the semi-batch column 
runs using toluene volume such that nearly all was 
caught the turbulent dispersion process. here successive 
data are listed Table for these runs, which lasted either 
five ten minutes, these apply the start run, the 5-minute 
mark, and the end run. The discrepancies successive 
values for given run can accounted for the basis 
experimental equilibrium data and 
the finite time required develop concentration boundary layers 
the non-flowing toluene phase. The last factor would explain 
the consistently downward trend the efficiency values with 
time near the start run. The two three values for 
given run were given equal weight including them all 
subsequent plots. 

Comparison carboy run with column run both using 
spout pipe i.d. 0.180 in. and zero toluene flow, indicates 
significantly higher stage efficiency the column average) 
than the carboy (45%). may argued that much all 
this difference can accounted for the larger water 
flow rate the column, where the water rate was 0.6 
gal./min. opposed 0.5 U.S. gal./min. the carboy. 
carboy run with column runs and all for 

spout pipe i.d. 0.269 in. and zero toluene flow, indicates 
however that even for the same water rate 1.0 gal./min., 
the column efficiency (65% average) considerably greater 
than the carboy efficiency (42%). occurs, moreover, 
despite the fact that spout heights were lower the column 
than the carboy. The principal explanation for this difference 
probably the fact that the average residence time droplet 
the dispersion zone the considerably greater 
than that droplet falling the toluene-water interface 
the carboy. column runs and above, the height the 
dispersion zone was in., while carboy run the spout 
height was in. For drops similar size falling equal 
velocities through similar medium, simple mass_ transfer 
dictates that for the settling process, 
directly proportional the height. this case the 
for the overall carboy process approximately 1.9, which 
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Figure 10—Effect water spouting rate stage efficiency 
static toluene, using schedule-40 spout pipe 
4-in. column. 


not much greater than the settling height ratio This 
near agreement gives support the suggestion that droplet 
residence time important the overall mass transfer process. 
addition, the added turbulence the dispersion zone the 
column probably increases the overall mass transfer coefficient 
around the droplets, while the stronger downflow currents 
toluene probably increase the mass transfer coefficient about 
the rising spout, over that the carboy. 

Though the settling droplets probably contributes sub- 
stantially the overall mass transfer process, there evidence 
that not the sole contributor. This evidence provided 
Figure which shows the effect toluene flow rate 
extraction efficiency for two constant water rates. droplet 
settling were all-important, the effect increasing toluene rate 
would increase mass transfer the water drops, which 
are moving net direction opposite that the toluene. 
fact, however, there slight decrease Eyg with toluene 
rate, particularly the lower water rate. This can only 
explained the contention that significant portion, though 
not necessarily most, the mass transfer occurs during the 
ascent the spout and fountain rather than during the descent 
the droplets. noted previously the discussion 
Figure increasing the toluene rate decreases the net relative 
velocity between the water spout and the cocurrent toluene, 
especially the lower water rate. This causes decrease 
overall mass transfer coefficient between the water spout, the 
height which increases only slightly (see Figure 4), and the 
cocurrent toluene, and hence apparently small decrease 
extraction efficiency for the overall process. 

The effect water spouting rate Eyg, using 1/8-in. 
schedule spout pipe and toluene flow, shown 
Figure 10. Increasing the water rate obviously has most 
decided positive effect Eyg about 1.2 US. gal./min. 
This can explained 
(1) the increasing spout velocity which increases the overall 

mass transfer coefficient about the rising spout; 
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Figure 11—Effect spout pipe diameter stage effi- 
ciency static toluene, operating near overhead entrain- 
ment point 4.in. column. 


(2) the increasing spout height (Figure 3), which counteracts 
the lowering time per unit length spout; 

(3) the increasing formation droplets 5), giving rise 
more mass transfer surface; and 

(4) the increasing depth dispersion zone with water rate, 
trend observed also Ryon their settlers, 
which results raising the residence time the droplets. 

Above 1.2 gal. water/min., factors and some extent 

taper off, with the result that tapers off sharply. 

gal./min., overhead entrainment water droplets starts 

(Figure 3). The plateau between 1.2 and 1.5 gal./min. 

obviously the desirable operating region, using standard 

spout pipe. 

The effect spout pipe diameter operating at, 
within 20% of, the overhead entrainment point each case, 
presumably within similar flat plateau region, shown 
Figure 11. The initial increase with spout pipe inside 
diameter undoubtedly due the very large initial increase 
with spout pipe diameter the allowable water rate before 
overhead entrainment will occur. This large increase volu- 
metric water rate more than counteracts several opposing 
tendencies, which only come into play decisively once 
increase spout pipe diameter longer accompanied 
large increase allowable water rate. These tendencies, 
which result the subsequent decrease with spout pipe 
inside diameter, are the diminution spout inlet velocity, 
spout height, droplet formation and dispersion zone depth 
with increase spout pipe diameter. Clearly optimum 
spout pipe inside diameter for mass transfer about 0.25 in. 
indicated Figure 11. 

Capacity Limitations 

picking optimum operating spout pipe diameter, 

have not only consider mass transfer, but must also take 


into account capacity limitations. Capacity restricted 
operating instabilities caused entrainment, loading, and 


OVERHEAD ENTRAINMENT OF WATER IN ne 


MINIMUM WATER RATE FOR ENTRAINMENT, U.S. GPM. 


o2 Q3 04 os 


a7 08 


SPOUT PIPE 1.0, in 


Figure 12—Initial entrainment points for water spouting 
static toluene (overhead entrainment requires slight 
toluene flow). 


eventual flooding one phase the other. Initial load points 
are shown Figure 12, which based largely data from 
Figure spout pipe diameter increases, the water rate 
necessary produce droplets sufficiently small entrained 
overhead slight toluene flow increases. the other 
hand, bottom entrainment toluene water independent 
spout pipe diameter, and, under the severe conditions imposed 
the contracting neck the carboy, occurs water rate 
low 1.2 gal./min. such unnecessarily severe 
conditions are imposed, then the maximum recommended spout 
pipe diameter, the intersection the two entrainment lines, 
would about 0.23 in., value which, sheer coincidence, 
very close the 0.25 in. obtained optimizing mass transfer. 

Actually the column, due the transition piece the 
bottom, the outlet water was slowed down that 
bottom entrainment toluene did not occur below about gal. 
water/min. above this transition piece, however, there 
was tendency for bottom entrainment occur for water 
rates considerably excess 1.2 gal./min., the special figure 
for the carboy. Thus, which involved neither end 
contractions expansions, spout pipe inside diameter 
0.23 in. could safely exceeded, certainly the optimum 
value 0.25 in., even 0.269 in., standard 1/8-in. 
schedule pipe. 


design for spouting water into toluene continuously 
were based 1/8-in. schedule spout pipe, water 
flow 1.2 U.S. gal./min., which the beginning the 
plateau Figure 10, then the corresponding toluene rate should 
such that any possibility flooding precluded. Using the 
flooding correlation for spray towers Minard and 
simplified applying the generous safety factors 
recommended, toluene rate not exceeding U.S. gal./(min.) 
(sq. ft.), corresponding 35% the calculated rate for water 
rejection, indicated. the 4-in. column this figure represents 
toluene flow 1.6 gal./min. 


Design Proposals 


One conception larger scale spouted mixer-settler 
illustrated Figure 13. The denser liquid introduced the 
bottom head through one more inlet ports and moves upward 
through several equally spaced spout pipes parallel, while 
the lighter liquid introduced the top head and moves 
downward through several pipes which surround the spout 
pipes symmetrically. The corresponding separated phases leave 
through ports immediately above the bottom sheet and below 
the top sheet, respectively. The heavy outlet phase passes over 
sliding weir, which used for interface control. 
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Figure Industrial-scale spouted mixer-settler. 

Toluene-phase inlet; Water-phase inlet; Toluene-phase 

outlet; Water-phase outlet; Adjustable weir; Bottom 

and top heads; Bottom and top sheets; Spider sup- 

port; Vertical baffles; Toluene-phase inlet pipes; 
Water-phase spout pipes. 


order gain the enhancement stage efficiency appar- 
ently contributed confining walls, vertical baffles are provided 
the primary design completely encompass the dispersion 
zone. These baffles partition off each spout from its neighbours 
and are mounted the outside walls the unit. typical 
compartment, shown Figure 13, would contain 
central spout pipe surrounded two four light phase inlet 
pipes, all held rings and spider supports from the baffles. 
Nozzle injectors could replace simple spout pipes provide 
further insurance high stage efficiency, but only greater 
initial and operating cost. 


assign feasible dimensions the above unit, the example 
water spouting into toluene considered. For this system 
has been shown previously that 1/8-in. schedule spout 
pipes, operating with water flow 1.2 U.S. gal./min. per 
spout pipe and toluene flow gal. /(min. (sq. ft.), 
are indicated. The last figure conjunction with the desired 
toluene throughput would determine the total cross-sectional 
area the extractor, while the water rate per spout pipe 
conjunction with the desired water capacity would determine 
the number spout pipes. The toluene inlet pipes not have 
critical size and could 1/4-in. schedule 40, while 
compartment, the light the incomplete data wall 
effect provided the present study, could reasonably four 
inches square. The distance between bottom and top sheets 
set ft., thus allowing ft. for the dispersion zone, 
though only ft. are anticipated, and least 1.5 ft. above and 
below this zone. The vertical baffles are ft. high, leaving one 
foot free space top and bottom permit removal the 
and extract phases. The water and toluene inlet pipes 
each extend 3.2 ft. from their respective sheets and therefore 


overlap almost five inches, which ample vertical length 
for the ring supports. 


should noted that quite possible and feasible 
reverse the direction spouting any all compartments, 
that is, have the present case downward spouting toluene 
into water. the above design, would necessary only 
lower the toluene inlet pipes compartment the same 
level the original spout pipe, and raise the spout pipe 
the same level the original toluene inlet pipes. 


Another possibility use array toluene 
and water inlet pipes, whereby the turbulence fostered the 
presence spout- -confining substituted the turbulence 
created the interaction adjacent spouts. This would 
eliminate partitioning baffles and cut down the necessary 
number toluene inlet pipes, hence effecting economy. 
illustrated Figure 13, the ratio the number toluene inlet 
pipes the number spout pipes can reduced from 4:1 
(or 2:1) 1:1 use unbaffled array. Whether the stage 
efficiency such array would match that baffled unit 
not, howev er, certain. 

The unbaffled array could also lend itself two-way spout- 
ing, water into toluene and toluene into water this case. 
increased vertical allowance, perhaps foot, would required 
for the dispersion zone. very likely that two-way 
spouting, which bears superficial resemblance Russian 
two-way injection could more than match the stage 
efficiency baffled one-way spouting. 


Comparisons and Conclusions 


The relative compactness extractor can measured 


and having the dimensions reciprocal time. Based the 
system water-benzoic acid-toluene and the recommended design 
conditions above, for which the fractional stage efficiency and 
therefore are estimated 0.90, for spouted mixer-settler 
works out This value, according Mathers 
pump-mix and their own air-lift contactor, 
and three times that obtained unpulsed packed 
all for the system water-acetic acid-methyl isobutyl ketone. 
Since the unpulsed packed column figure similar those 
shown the packed column data Appel and the 
system water-benzoic acid-toluene, the procedure comparing 
contactor effectiveness using data for the two different systems 
considered legitimate. 

cascade spouted mixer-settlers for multistage operation 
would require inter-pumping the phases. Making liberal 
allowances for various frictional head losses the present 
equipment, the energy requirement for the recommended desi sign 
litre/min. total flow. According Mathers and 
this value about one-quarter that displayed mechanically 
agitated but about four times that shown their 
own air-lift contactor. The degree mixing was comparable 
the three units, since the stage efficiencies were all close 
90% and the mixer retention time about 0.5 min. the spouted 
contactor the was taken the middle ft. dispersion 
zone. 


spouted mixer-settler has thus been shown 
tive continuous-flow unit, which operates without the large 
energy losses inherent mechanical agitation, performs both 
mixing and settling single relatively compact chamber, 
gives high stage efficiency when opt timally operated, and has 
very simple design. Although its energy consumption 
greater than that air-lift has the advantage 
over the latter not producing any air contaminated with 
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entrained liquid, advantage which pertinent 
chemical work. 

Liquid-liquid spouting may characterized bearing 
similar relationship conventional relatively bulky spray tower 
operation up-down prilling does the older bulkier one-way 
prilling, used for example the manufacture ammonium 
nitrate fertilizer granular 
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Nomenclature 

bulk concentration solute, ft. 

value for toluene phase, ft. 

value for aqueous phase, ft. 

Cr* value equilibrium with Cy, lb.-moles/cu. ft. 

fractional stage efficiency, defined Figure 

contactor effectiveness, defined Equation (4), 

Murphree extract efficiency, defined Figure 

Murphree spouted-phase efficiency the present runs. 

Murphree raffinate efficiency, defined Figure 

total flow rate both phases, cu. ft./hr. 

spout height, ft. 

number theoretical stages for given operation 

temperature, °F. 

space-average velocity spouted liquid leaving spout 
pipe, ft./hr. 

Umax Maximum local velocity spouted liquid leaving spout 
pipe, ft./hr. 

total volume contactor (mixers plus settlers), cu. ft. 

kinetic-energy correction factor take account radial 
velocity distribution spout pipe. 

given value 

density continuous-phase liquid, ft. 

density spouted liquid, ft. 

Subscripts 

entrance 

exit 

temperature 
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Mass Transfer from Fluid and Solid 
Spheres Low Reynolds 


and TRASS? 


Predicted continuous phase mass transfer co- 
efficients for fluid and solid spheres the Stokes 
region flow are correlated Sherwood numbers 
plotted against Peclet numbers. These correlations 
permit study the effect internal circulation 
the rate mass transfer from fluid spheres. 
Peclet number the theory predicts three- 
fold increase transfer due the circulation move- 
ments within gas bubble. Peclet number below 
10°, the difference between circulating and non- 
circulating (or solid) sphere indiscernible. 


The analysis presented this paper for fluid 
transfer around solid spheres. more elaborate 
concentration profile through the boundary layer has 
been used the present investigation and the result 
reported for solid spheres differs about from that 
given Friedlander. 


The more interesting details the mathematical 
treatment and the machine computation the result- 
ing differential equations are included. 


rate mass transfer from spherical particles moving 
uniformly continuous fluid medium has direct application 
many types two-phase contacting 
nately precise mathematical formulation this problem 
exceedingly difficult the flow regime practical significance. 
However some theoretical studies diffusion from solid spheres 
low Reynolds numbers have been carried 


This paper mathematical analysis the influence 
the internal circulation predicted the transfer 
rate from spherical fluid particles Specifically, the 
continuous phase mass transfer coefficient was calculated for 
typical liquid-liquid and gas-liquid system and compared with 
that obtained for the analogous solid sphere case. This work 
was suggested part the observation Hammerton and 
that the transfer rate from circulating bubbles 
slightly soluble pure gases was about five times large from 
non-circulating bubbles. This increase cannot attributed 
improved mixing inside the sphere (since the entire resistance 
transfer the continuous phase) and therefore must 
due the effect internal circulation the external flow 
received May 14; accepted October 10, 1960. 


2Department of Chemical Engineering, University of Toronto, Toronto, 
Ont. 
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Diffusion Model 


The model used this work similar that used 

assumed that steady state diffusion occurring from 
fixed spherical particle radius ‘a’ continuous fluid medium 
laminar flow. 

Two cases will considered both which have resistance 
transfer present the external fluid only. Case illustrated 
the left side Figure represents diffusion from solid 
sphere where the velocity distribution the continuum 
obtained from the solution the linearized Navier- 
Stokes equations. Case the right side Figure considers 
diffusion from two-component fluid system here the circulat- 
ing disperse particle saturated with the continuous phase and 
where the velocity distribution the continuum obtained 
from the Hadamard™ solution the linearized Navier-Stokes 
equations. 

concentration boundary layer defined having radius 
which function the angle with the concentration 
varying from the saturation value the surface the main 
stream value the outer edge. Natural convection neglected. 


Mathematical Relationships 
The amount solute passing any section indicated 


stream function, the volumetric rate flow, and 
the concentration with respect main stream concentration 
zero. 

Assuming molecular diffusion the direction negligible, 
the difference this quantity between any two sections separated 
the angle can equated the amount diffused from the 
surface expressed Fick’s first law. Thus 


dy’ 


When the non-dimensional symbols Friedlander are used: 


, 


a 


then Equation (1) becomes: 


Ww 


STOKES STREAMLINES — 
CONCENTRATION --- 
BOUNDARY LAYER 


SOLID SPHERE-Pe=100 


DIRECTION FLOW 


HADAMARD STREAMLINES — 


CONCENTRATION 
BOUNDARY LAYER 


GAS BUBBLE - Pe=100 


Figure 1—Comparison solid and gas sphere forced 
convection (Re 1). 


The important quantity this study the amount diffused 


from the total surface the sphere which given 


o 
180°. This integral, evaluated the methods discussed 
the next section, contains the unknown boundary layer 
thickness However, further substitution into Equation (3) 
gives rise ordinary non-linear differential equation which 
can solved numerically for 

solve this integral, first necessary express both 
(a) Stream Function (r, 

The linearized form the Navier-Stokes equations can 


expressed as: 


This equation has the general solution 


The specific solution for solid spheres obtained 
follows: 


non-dimensional form: 


The boundary conditions for circulating fluid sphere 
given were applied Equation (5) giving 
the following: 


where and are the continuous and disperse phase viscosities 


respectively. Since Equation (8) represents the stream function 
for particle moving quiescent liquid was necessary 


superimpose the flow field negative velocity equal the 
Hadamard velocity fall. The resulting equation written 
non-dimensional form is: 


For simplicity, both Equations (7) and (9) can written 


where for Equation (7); for Equation (9) 


for Equation (7); for Equation (9) 


(b) Concentration (r, 


mass balance elemental cube the boundary layer, 
considering bulk flow both the radial and tangential directions 
and diffusional flow the radial direction, gives rise the 
following equation polar coordinates. 


Oc’ 
= 


where and are the radial and tangential 
spectively. 

This partial differential equation leads required 
but the solution quite complex. obtain expression 
more suitable for the subsequent calculations, arbitrary 
concentration polynomial was adopted which satisfied known 
boundary conditions. These conditions were follows: 


the above, (1), (3) and (4) came directly from the 
physical model. Condition (2) was obtained applying 
Equation (11) the sphere’s surface. will noted that since 


atr 


atr 
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and are zero the surface for the solid and fluid 
spheres considered this work, Equation (11) reduces (2) 
for both cases. 

The various polynomials investigated are shown Table 
Most did not satisfy one more the following additional 


(a) boundary conditions must independently satisfied. 
(b) discontinuities can occur within the boundary layer. 


corresponds diffusion into stagnant fluid and must 
satisfied. 


noted column four Table only polynomials (2), 
(5) and (10) met these added conditions. 


Polynomial (10) satisfied the four boundary conditions and 
was therefore used the subsequent calculations. addition 
polynomial (2) was also carried through final solution 
two-term profile was sufficiently accurate characterize 
the concentration the boundary layer. This latter profile 
the one used Friedlander (in slightly different form) for 
mass transfer from solid spheres and his work served guide 
and comparison for that portion the present study. 

The concentration profiles obtained from polynomials (2) 
and (10) are respectively: 


2r(y 


(c) Evaluation 
simplify the boundary limits, Friedlander transformed 


TABLE 
Boundary 
No. Form Polynomial Conditions 
(1), (3) (c) 
c=a (1), (3) 
(1), (3), (4) 
(1), (2), (3) (b) and (c) 
(1), (3), (4) (c) 
(1), (2), (3) (a) 
| 


Boundary conditions must independently satisfied. 
(b) discontinuities can occur within the boundary layer. 


(c) 


TABLE 


Profile Profile 


(ey), 


Substituting for the form Equation (10) and from 
Equations (13) and (14), the following two expressions are 
obtained: 


Using Equation (13) 


Using Equation (14) 


For convenience later computations, these equations were 
simplified 


The coefficients were evaluated for the three cases 
solid, liquid and gaseous spherical particle and are given 
Table The specific examples were chosen be: 


(a) solid (in any fluid) 
(b) liquid acetate drop rising water) 
(c) gas (carbon dioxide bubble rising water) 


Determination Mass Transfer Coefficient 


Equations (18) and (19) when substituted into Equation (3) 
give respectively the following differential equations general 
form: 
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| 


Now Fick’s first law written dimensionless form becomes: 


and substituting from Equation (3), one obtains 


The final results will reported this form. 


(a) Analytical Solutions Equations (20) and (21). 


Friedlander his sphere analysis (using profile showed 
that certain order magnitude approximations were possible 
for thick and for thin boundary layers and, for these limiting 
cases, obtained analytical solutions Equation (20). The same 
techniques were applied the present study Equation (21)— 
again the solid sphere case. These analytical solutions are 
shown Table general analytical thin boundary layer 
solutions were possible for fluid spheres. However the thick 
boundary layer solutions for both solid and fluid spheres reduced 
This would expected since very thick boundary 
layer implies diffusion into stagnant medium. 


(b) Numerical Solutions Equations (20) and (21). 


Except the restricted regions mentioned above, these 
equations required numerical solution. The method used 
the authors discussed below some detail several interest- 
ing points arose. The entire calculation was programmed 
the 650 computer. 


Both Equation (20) and (21) can expressed the follow- 
ing form 


where and are given Table 


TABLE 


Profile Profile 


Thin 


Pe) 


Solid and Pe+— Pe?+... 

Pe) 


From physical considerations dy/d@ must positive for the 
infinite unless (2/Pe cos also vanishes which 
case indeterminate. Thus the solution the equation 
passing through the point will follow the line 
reached which solves the equation 


this point, the slope the required solution will change 
discontinuously and will swing into the region Thus 
solution cannot initiated. The first step was therefore 
solve Equation (25) for (the Newton-Raphson iterative 
was found convenient for machine computation) and 
start the solution from the point y*, 

The numerical computation itself used power series 
(Equation (27)) obtain four starting values (at 0.5, 
1.5, 2°) and employed these 
formulae suggested The frequently used 
Milne method proved unstable for this problem. 

For convenience, the substitution cos was made 
Equation (21) giving 


The solution was assumed have the form 


Also and can expanded Taylor Series follows: 


92”)? 


Substituting Equations (27), (28) and (29) into Equation (26) 
and collecting like powers one obtains the following 
expressions for and 


9 


After the four starting values were obtained, predictor 

Equation (31) and corrector Equation (32) given below were 


applied. The corrector equation was repeated until two suc- 
cessive values agreed four significant figures. 


The increment was the angle radians corresponding 4°. 


a= + fo’) ~~ 


TABLE 


Ce 


Using Equation (13) (Profile 


Using Equation (14) (Profile 10) 
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|. CIRCULATING GAS BUBBLE- CARBON DIOXIDE IN WATER 
2. CIRCULATING LIQUID DROP-ETHYL ACETATE IN WATER 
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RESULTS USING PROFILE 10 
———— RESULTS USING PROFILE 2 


Sherwood Number 


Peciet Number 


Figure 2—Transfer from fluid and solid spheres. 


this manner values were obtained intervals 
tion, the was also obtained (using Equation 19) and 
the value the integral 180° estimated extrapolation. 
This latter quantity was then substituted into Equation (23) 
obtain the Sherwood number and hence the mass transfer coefh- 
cient. The results for solid, liquid and gas particles for both 
concentration profiles are shown Figure 


Discussion and Conclusions 


The results show that Peclet number circulation 
movements within gas bubble increase the mass transfer rate 
about threefold over that non-circulating bubble solid 
particle. The enhancement typical circulating liquid drop 
the same Peclet number about two and one-half times that 
solid particle. For any particular system, the exact degree 
improvement depends the viscosity ratio the disperse 
and continuous phases. The increase due circulation decreases 
with decreasing Peclet number and almost vanishes for Pe< 
This expected since low Peclet number high 
diffusivity) indicates very thick concentration boundary layer 
where the relative effects circulation are much 
nounced. 


Hammerton and reported that fivefold reduction 
mass transfer occurred when ethylene bubbles rising water 
had previously contacted petroleum contaminant 
presumably prevented circulation. Since their work was 
much higher Reynolds number than that considered the present 
study, quantitative comparisons are not possible. However the 
general agreement with the results reported above suggests that 
circulation higher Reynolds numbers may play similar role. 
should pointed out that circulating fluid sphere has 
higher Peclet number due its increased terminal velocity and 
hence some increase in transfer rate over its non- -circulating 
partner would expected. The curves the present investiga- 
tion, however, are independent this effect since they have 
been plotted against the Peclet number. 


Part the present study was evaluation various 
assumed concentration profiles and their effect the final 
solution. shown Table many the polynomial forms 
attempted gave profiles which were consistent with the physical 
picture. The four term profile (No. 10) slightly different 
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Sherwood-Peclet number relationship than the two term profile 
(No. 2). The greater complexity involved the former must 
weighed against the accuracy desired, taking into account the 
inherent limitations the defining model. effects such 
tangential diffusion, natural convection and interfacial tension 
could mathematically included, the more complete profile 
would undoubtedly warranted. solid sphere 
heat transfer study for Peclet numbers less than obtained 
theoretical results identical those given here for profile 10. 


The quantity y*, although its use here was overcome 
mathematical hurdle, can given the physical significance 
back diffusion from the forward stagnation point into the 
on- coming fluid stream. Since infinite concentration boundary 
layer realized the number zero, this hypo- 
thesis seems reasonable. Indeed approaches infinity for this 
condition can seen from Equation (25). 

Further details the mathematics derivations given this 
paper may found either reference (10) (11). 


Nomenclature 

sphere radius, cm. 

constants concentration polynomial 

A,B constants stream functions 

dimensionless 

concentration, g.-moles/cc. 

saturation concentration, g.-moles/cc. 

diffusivity, cm.?/sec. 

mass transfer coefficient, cm./sec. 

Peclet number, dimensionless 

dimensionless 

Reynolds number, dimensionless 

Schmidt number, dimensionless 

Sherwood number, 2ka/D, dimensionless 

main stream velocity, cm./sec. 

radial velocity, cm./sec. 

tangential velocity, cm./sec. 

distance direction flow, cm. 

dimensionless 

boundary layer, cm. 

angle from negative x-axis 

continuous phase fluid viscosity, gm./cm.-sec. 

disperse phase fluid viscosity, gm./cm.-sec. 

stream function, dimensionless 


stream function outer edge boundary layer, 
dimensionless 


stream function, cc./sec. 
constants 
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Dynamic Analysis Bubble Plate 


frequency response method examined 
possible experimental method for the simultaneous 
determination liquid mixing, plate efficiency, 
vapor-liquid equilibrium, and liquid holdup 
single bubble tray. The method depends upon the 
introduction sinusoidal fluctuations the concen- 
tration species the liquid and vapor fed the 
plate and the measurement the induced concentra- 
tion variations the same species the effluent 
liquid. applicable multi-component systems. 


The analytical results used interpret the 
measurements are based linearized model. 
Despite this idealization, systems exhibiting marked 
the liquid-vapor equilibrium relation 
are adequately represented. Systems showing very 
significant enthalpy effects are treated quantitatively 
when the liquid the plate perfectly mixed. This 
done with the aid the “signal flow graph” 
technique. 


With incompletely mixed plate liquid, mathe- 
matical difficulties have far prevented the quan- 
titative interpretation measurements the fre- 
quency response when enthalpy effects are large. 
present appears that the frequency response tech- 
nique limited perfectly mixed plates, i.e., 
small laboratory sized plates, incompletely 
mixed systems which the heat effects are small, 
i.e., gas absorption plates and distillation systems 
which the volatilities and latent heats all com- 
ponents are comparable magnitude. 


advances the fields process design and 
control have increased the need for more complete infor- 
mation concerning the behavior chemical engineering systems. 
particular, now recognized that the optimization the 
design and operation chemical plant requires the specifica- 
tion its transient behavior. The information needed predict 
the unsteady-state performance plant component not 
easy obtain and faster, cheaper methods than now available 
are much desired. The motivation the work reported 
here stemmed from the need predict the dynamic response 
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bubble plate contactors, class equipment which found 
all branches the chemical process industry. 

this paper, the feasibility using frequency response 
technique the simultaneous experimental determination 
liquid mixing, plate efficiency, vapor-liquid equilibrium and 
liquid holdup single experimental bubble tray investigated. 
The method involves the introduction sinusoidal composition 
fluctuations the concentration one the species entering 
the plate with the liquid and vapor feed and the measurement 


the composition fluctuations that same species the 


liquid from the plate. The method applicable multi-com- 
ponent systems. 

Well developed techniques exist for the analysis the 
dynamic response linear systems and the form the perturba- 
tion imposed linear system for the purpose determining 
its dynamic behavior relatively unimportant. Neither 
these circumstances apply the case non-linear system 
and unfortunately bubble plate contactor non-linear device. 
For small fluctuations, however, linearized model can 
expected approximate the behavior the actual plate. Under 
such circumstances, the correspondence between the behavior 
the actual plate and that predicted the basis linearized 
model needs known function the amplitude and 
type disturbance imposed. the present paper, the form 
the disturbance limited sinusoid and attention focussed 
upon the information that can obtained from applying such 
disturbance actual plate. 

Two important considerations led the choice sinu- 
soidal disturbance: 

Should prove capable producing useful results, 
sinusoidal forcing technique would important experimental 
tool. Firms concerned with the design and control separation 
equipment using bubble contacting plates might find profitable 
build both laboratory and commercial sized plates which 
critical separations could studied means the measure- 
ment the response the plate sinusoidal oscillations 
the input streams. 

operating plate can ascertained from measurements 
the response the plate the fluctuations input 
quantities which occur normal part its commercial 
operation. this were possible, artificial disturbances 
would required. However, present the feasibility 
such technique rests the use correlation techniques 
which require that the response the system sinusoidal 
perturbations understood. This suggests that the response 
the system sinusoid variations studied before investi- 
gation the utility random fluctuation technique under- 
taken. 
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The Linearized Model 


Consider single bubble plate. Idealize the plate follows. 


Represent the liquid mixing the plate the Kirsch- 
baum idealization series perfectly mixed 


each pool moles liquid shown Figure 


The liquid flows from one pool the next without bypassing 
pools back-mixing. Although this simplification the mixing 
process gross approximation the actual state affairs, 
shown Reference that the range customary com- 
mercial operation, the pool concept essentially equivalent 
such alternate mixing hypotheses 
and 


The vapor enters the plate with uniform composition 
and divides equally among the pools. 


The liquid flow rate and the vapor flow rate V/N are 
unchanged passage through pool and are invariant with 
time. important note that this assumption ignores the 
fact that unequal heats vaporization and sensible heat effects 
can change the vapor flow rate the expense the liquid flow 
rate. Later, will shown that such energy effects can 
significant and may have important influence the dynamics 


. 
plate behavior. 


assumed that the vapor composition equilibrium 
with the liquid can represented linear relation over the 
composition range found the plate: 


not assumed here that the vapor-liquid equilibrium relation 
truly linear. Rather, that over limited composition range, 
the true curve can approximated its cord. and then 
are actually functions composition which can approximated 
average values for small variations Later, the validity 
this approximation explored. 


The composition the vapor leaving pool ys, 


and assumed constant across the plate. 


With the above idealizations, mass balance written about 
pool for the same component gas and liquid gives 


Equation (3), the subscript refers the pool. and 
refer chemically identical species the liquid and vapor 
respectively. Equations the form and apply each 
component the system and restriction placed the 
number components present. 


eliminating and from Equations (1), (2), and (3) 
there results 


ri 
2\<! 
2\<! 


2\<) 


Figure 1—The representation liquid mixing series 
perfectly mixed pools. 


mEV dx, 
When the inlet vapor and/or liquid composition 
the sum mean value and time-dependent value 
steady-state operation the plate. the above substitu- 
tions Equation (4) becomes 


mEV 


the steady-state expression, 


mEV 


subtracted from Equation (5), the relationship amongst the 
time-dependent values found 


Taking initial conditions Equation (1) steady-state 
operation, then the initial conditions Equation (6) are 
x’; = x’) = x's x'y = = (7) 


The Laplace transform Equation (6) with the initial 
conditions(7) 


Equation (8) linear finite-difference equation with constant 
whose solution readily found 


VEm 
(Hs VEm) VEm 


Equation (9) the desired result. When replaced 
and jw, the result expresses the response the outlet 
liquid composition sinusoidal composition variations 
either the inlet liquid the inlet vapor 


used denote the ratio the amplitude 
the outlet liquid composition variation the amplitude 
sinusoidal variation frequency imposed the entering 
liquid composition when the inlet vapor composition and all 
flow rates are maintained constant, Equation (9) gives 


Rearrangement Equation (10) yields 


This statement predicts that several frequency response 
experiments are performed and the results plotted 


2 


vs. straight line formed whose slope and 
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H dx’, 5 

N N | 


mEV 

must plotted for several values and that one which 
yields the best straight line selected the correct one. With 


—\2 


determined. There remains only separate and from 
the product 

Paranthetically, should mentioned here that Equation 
(11) implies method determining liquid holdup and liquid 
mixing imperfectly mixed bubble plate sinusoidal 
variation the concentration nonvolatile solute 
easily seen that although Equation (11) was derived for 
volatile component, can applied equally well nonvolatile 


Since not known, curves 


one simply ignoring the term mEV the expression for the 
intercept. The nonvolatile may thought 
component whose vapor-liquid equilibrium curve horizontal 
line coincident with the axis. The liquid holdup 
and mixing parameter are determined described above, 
guessing plotting vs. and finding the slope 


Return now the case volatile components and the 
problem separating and from the product mE. Let 
represent the ratio the amplitude the outlet liquid 
composition variation the amplitude sinusoidal variation 
frequency imposed the entering vapor composition 
when the entering liquid composition and all flow rates are 
held constant. Then (9) gives 


where 


Three alternative methods may used obtain the values 
from the procedure represented Equation (12). 


First, determined sufficiently high frequency 
such that 


then Equation (12) implies 


With knowledge and from the results set 
high range sufficient determine from Equation (14). 

possible obtain more than one value 
different frequencies, then independent estimates both 
and are available plotting vs. since the inverted 
square Equation (14) 


straight line fitted such plot should have intercept 


and slope from which both and could 


calculated. 
determined the same mean conditions and 


shows that Hence, the frequency range required 


the output response liquid input composition variation. 

The disadvantage the above method lies the fact that 
high frequencies, well becomes small 
and the actual measurement increases when 
variation must detected precisely. relatively 
simple experimental techniques are available for detecting the 
amplitude and phase sinusoid known frequency the 
presence quite high level noise and distortion. 

The second method extracting information from Equa- 
tion (12) involves experiments with very low frequencies long 
sinusoidal periods. seen that when substituted into 
Equation (12) and evaluated vanishingly small values 


there results 


This implies that curve vs. extrapolated 
yields the value 


then known the square root the intercept 


Equation (11), value can determined. 
third approach available neither sufficiently high nor 


low frequency experiments are practicable. This involves 
direct determination the absolute value Equation 
quency, plot 


yields straight line, 


from which and are easily found. order find 


both the absolute value and the argument the complex 
number must known functions frequency. The 
absolute value simply The argument may found 
Moivre’s theorem. For 


= 


and 


where 


With this, Equation (17) becomes 


V 
The values and needed for the evaluation 
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the angle must obtained from determination the 
slope and intercept Equation (11). 

The validity Equation (18) depends upon the requirement 
that Equation (17) has the same values the which 
replaces it, the transfer function liquid outlet composition 
with respect liquid inlet composition. This means that the 
two amplitude ratios which are combined the calculation 
the left-hand member Equation (18) must the result 
pair experiments which the 


mEV 


equal for both. This implies that the frequency and mean 
value all input variables both experiments should the 
same. 


Testing the Linear Model 


The previous discussion makes clear that the linear 
model obeyed the frequency response method indeed 
powerful experimental tool. 

preliminary evaluation the linear model was carried 
out comparing the linear model predictions with the results 
digital computer simulation operating bubble plate. 

Computer simulation not replacement for experimental 
test using operating plate. Rather, very useful means 
initial evaluation and exploration. Wide ranges operating 
and physical parameters may covered under conditions that 
facilitate interpretation their effect. the linearized model 
shows satisfactory correspondence with realistic computer 
simulation actual plate, then necessary and desirable 
test the method means physical measurements 
operating plate whose characteristics have been determined 
conventional procedures. 


Non-Linear Vapor-Liquid Equilibria 


One the most questionable assumptions required the 
linear model that linear vapor- liquid equilibrium relation- 
ship over the range the plate composition. Consequently, 
simulation was built which included unusually severe local 
curvature the vapor-liquid equilibrium line. The equilibrium 
The quadratic 


line chosen was 0.60 Figure 
expression 


VAPOR MOL. FRACTION, 


0.5 1.0 


LIQUID MOL. 


Figure equilibrium curve for Series 
and 


coincides with the equilibrium line the solid segment the 
sharp bend near 0.55. equilibrium line this form 
corresponds incipient phase separation. The simulation was 
run the IBM 704 computer the Computation 
Center and was subjected sinusoidal composition disturbances 
required the linear analysis methods. The digital response 
waves obtained from the simulation were analyzed harmonic 
content the discrete Fourier approximation method described 
The details the simulation program and 
complete numerical results are available Reference 18. 

Inlet conditions the plate were chosen that the liquid 
composition range across the plate also coincided with the bend 
Figure The calculations were carried out using 0.60, 
The steady-state liquid composition for each pool was computed 
using the actual equilibrium curve (Equation (19); the results 
are given Table 

Two sets computer simulation sinusoidal forcing 
the input streams were carried out. Series employed 
constant excitation frequency 0.632 radians/second and 
total six calculations were made. Three involved variation 
the inlet liquid composition 0.5%, 5.0% and 20.0% 
normalized amplitude about 0.60 with held constant 
0.70. Three involved fluctuation the inlet vapor 
composition 0.5%, 5.0% and 20.0% normalized amplitude 
about 0.70 with held constant 0.60. The results 
are given Table 

The simulation Series was made varying frequency 
was varied from 0.2 3.0 radians per second. ‘Six calculations 
were made employing fluctuation about 0.60 with 
held constant 0.70 and six calculations were made employ- 
ing fluctuation about 0.70 with held constant 
0.60. The results are given Table 

Examination Table discloses two interesting results. 
The non-linear vapor-liquid equilibria generates mean value 
the reponse wave which differs from the steady-state value. 
The differences are reflected the column marked 
harmonic amplitude”. The entries that column are the 
quantities Xoo the mean value 
the liquid composition leaving the plate (the 20th pool) and 
0.495725 the corresponding liquid composition leaving the 
plate under steady-state operation shown Table The 
harmonic analyses the results shown Table show that 
despite the highly non-linear equilibrium curve the wave forms 
are slightly distorted fundamental sinusoids. all cases the 
distortion chiefly due the second harmonic and the ampli- 
tudes the third harmonics are small comparison with those 
the first and second. result, the amplitude the 
fundamental wave, which the value required the linear 
analysis methods, can obtained without resorting Fourier 
analysis. good approximation may taken one- 
half the peak-to-peak distance the distorted wave. 


TABLE 


0.70, 20, EQuILIBRIUM ACCORDING EQUATION (19). 
Liquid Mol Liquid Mol 
Pool Pool Fraction 

0.587217 0.528004 
0.580804 0.523111 
0.568101 0.514060 
0.561878 0.509908 
0.555780 0.506003 
0.549834 0.502341 
0.544065 0.498918 
0.495725 
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TABLE 


Series AMPLITUDES THE HARMONIC COMPONENTS THE RESPONSE WAVES 


CONSTANT EXCITATION FREQUENCY 0.632 


Amplitudes Harmonic Components* 

nput 

No. Input 

Zeroth First Second Third Fourth 
0.050 0.002425 0.012226 0.000290 0.000003 
0.200 0.028356 0.061507 0.004946 0.000065 
0.200 0.001479 0.023770 0.000070 0.000001 


*These amplitudes are those the response wave 


(twentieth) pool. 


Figure the data Series Table are plotted 
accordance with Equation (11) for various assumed values 
the number pools, The results least squares 
evaluation the the results the straight lines predicted 
Equation (11) are given Table appears that the correct 
value can determined least squares treatment 
Equation (11) despite the non-linearity the equilibrium curve. 
Near the correct the variances are small and conceivably 
could masked experimental errors actual application 
but will shown this does not appear serious practical 


limitation. accurate value the holdup can found 
from the slope the correct line Equation (11) using 
results. The wrong lines shift the determination 
moderately. 


The results Table were used separate 
and All the methods discussed earlier were tested and 
were found give similar results. Only the phase angle tech- 
nique employing Equation (18) will covered here. Values 
were computed function using Table and several 
assumed values Figure shows the values plotted 
versus accordance with Equation (18). graph reason- 
able size cannot resolve the differences between the (on 
least squares basis) straight lines corresponding the differing 
assumed values and hence only one line shown. However, 
least squares analysis did give somewhat different values 
slope and intercept for different values and these results 
were used obtaining the values reported Table The 
columns Table are: (1) the assumed (2) the value 
obtained from the intercept the data treated 
Equation (18); (3) the liquid mol fraction calculated from 


x20 0.495725 


-—, where 0.495725 the steady-state mol fraction from the last 


Equation (19); (4) the point efficiency obtained from the 
—\2 


slope the data treated Equation (18) and using 


values from Table (5) transcribed from the 
column Table (6) the slope derived from and 


145 
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Figure 3—Series Least-squares straight lines 
Equation (11). 


TABLE 


MEAN VALUES THE RESPONSE AND AMPLITUDE RATIOS 


Normalized Input 


Run No. Excited Input Amplitude 
0.10 
0.10 
0.02 
D4 0.02 
0.02 
0.02 
0.20 
0.20 
0.06 
0.06 
Dit 0.06 
0.06 


Amplitude Ratio 


Frequency Mean Liquid 

3.0 0.495889 
2.0 0.494928 

1.0 0.495788 0.059970 

0.6 0.495836 0.217713 

0.4 0.338157 

0.2 

3.0 0.495752 

2.0 0.495877 0.034840 

1.0 0.495749 0.072533 

0.6 0.495800 

0.4 0.495911 0.187171 

0.2 0.496107 0.348895 
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TABLE 


2 


1.084950 0.350466 

1.079777 0.292278 

1.074580 0.250079 

1.069732 0.218187 

1.065289 
TABLE 


3.183 0.451 0.66358 0.83218 1.2541 0.531 


3.374 0.442 0.68193 0.78246 1.1474 0.536 
3.241 0.448 0.70056 0.73240 1.0455 0.540 
3.054 0.456 0.71950 0.68558 0.9528 0.544 
2.804 0.467 0.73856 0.64256 0.8700 0.547 


transcribed from Table 


seen from Table that satisfactory values the point 
efficiency are obtained even the value chosen differs 
from the correct one. the other hand, when the intercept 
Comparison the value Table with the steady-state 
values the liquid composition shows that the values fall 
below the range the plate. the other hand, the slopes 
calculated from the values obtained using the 
results and Equation (11) together with values from 
results correspond liquid mol fractions near the center 
the plate. 

From these results appears that even when the vapor- 
liquid equilibrium line shows marked curvature, frequency 
response measurements may coupled with analytical pre- 
dictions based upon linear model order obtain plate 
parameters the following procedures are used: and the 
product are obtained from experimental measurements 
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Figure 4—Series Least-squares straight lines Equation 
(18) for 16, 18, 20, 22, and 24. 


Sum Squares 


0.83218 14.2080 278 


eof vs. conjunction with Equation (11). The plat 
measurements used conjunction with either Equations (15) 
(18). The slope, the equilibrium curve computed 
from the value obtained from data alone and the 
average liquid composition the plate. 

The data coupled with Equations (15) (18) will 
also yield value this differs significantly from that 
found above, the curvature the vapor-liquid equilibrium line 
pronounced. 


The physical basis for the superiority the linear approxi- 
mation analysis using data due the following. 
When fluctuation introduced into the inlet liquid, the induced 
fluctuations must across the entire plate before reaching 
the outlet where they are measured. Consequently, the outlet 
fluctuations should reflect average parameters across the plate 
and particular, average slope corresponding average 
central composition. the other hand, when fluctuations 
are introduced into the incoming vapor order make 
measurements, the uniform distribution vapor across the 
entire plate causes the composition fluctuations introduced 


(B) 


Figure 5—Flow graphs and 
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(a) 


uniformly across the plate also. result, the effect the 
forcing fluctuation maximum the down-stream end 
the plate and the attenuation measured the outlet heavily 
weighted the characteristics the downstream liquid. 


The linearized frequency response method gives satisfactory 


values and even when the equilibrium curve exhibits 
marked curvature. The vapor-liquid equilibrium values obtained 
under such circumstances are less satisfactory. This primarily 
due incomplete mixing the plate and the resulting plate 
concentration gradient. well-mixed plate, the vapor- 
liquid equilibrium results approach the precision the 
and results. However, even with imperfect mixing the fre- 
quency response results may used estimate the vapor-liquid 
equilibria follows. Presumably experimenter does not 
know beforehand that the curve non-linear nor does know 
the liquid composition distribution the plate. However, the 
mean inlet and outlet liquid compositions, the mean inlet vapor 
composition, and the liquid and vapor flow rates are measured. 
From these parameters, steady-state mass balance around the 
plate determines the average effluent vapor composition, 
where average refers both time and plate area. Now the 
point efficiency defined 


However, the experimental plate the distribution and com- 
position the incoming vapor stream uniform and 
the point efficiency assumed uniform over the tray area. 
Hence, when Equation (2) averaged over both time and tray 
area the result 


still the point efficiency and determined from the dynamic 
analysis. Since are known, may calculated and 
provides the vapor-liquid equilibrium curve. 


The abscissa not known precisely there significant 
concentration gradient across the plate. single time mean 
composition sample taken just downstream the inlet sec- 
tion, comparison this composition with the outlet composi- 
tion will indicate the magnitude the composition gradient. 
inconvenient traverse the plate find the steady-state 
liquid composition distribution, the average the entrance 
sample composition and the outlet composition may serve 
reasonable approximation the area average liquid composi- 
tion. This composition the abscissa which sought for the 
point the vapor-liquid equilibrium curve. plotting curve 
from series such determinations, the values also 
determined the course the frequency analysis experiment 
fix the local slope the equilibrium curve and facilitate the 
construction the curve. 


Significant Latent and Sensible Heat Effects 


addition linear approximation the equilibrium 
curve, the linearized analysis has assumed that the liquid flow 


rate and the vapor flow rate are unchanged passage 
through pool. These flow rates can, course, change 
practice result energy effects. order ascertain the 
importance such effects, extreme case was examined. 
commercial scale extractive distillation plate separating 
hydrocarbons with the aid furfural-water stream was 
simulated the IBM 704 Digital Computer. The simulation 
included the actual vapor-liquid equilibrium curve and the 
and sensible heat effects associated with all streams. However, 
order understand the complicated interactions, the plate 


With the linearized relations Equations (10 and 
(12) become 


apparent that perfectly mixed plate for which the 
linearizing assumptions are valid not necessary run 
frequency response experiments over range sinusoidal 
frequencies order determine plate efficiency and vapor- 
liquid equilibria. single pair experiments the same 
frequency sufficient provided liquid feed composition 
fluctuated one experiment and feed vapor composition 
fluctuated the other, but the time-mean values all variables 
are the same both. Under these conditions, the ratio 
Equation (22) Equation (21) gives 


When the predictions Equations (21), (22) and (23) 
were compared with the performance the simulated extractive 
distillation plate, the correspondence was poor for and 
considered separately, but their ratio, Equation (23) 
was good agreement. order define the cause the 
discrepancies, theoretical analysis was made the plate 
behavior, taking into account the energy effects. The details 
the analysis are given the appendix. 

The principal finding that Equations (21), (22) and (23) 
should replaced 


Equations (24), (25) and (26) are approximations that certain 
the energy terms have been discarded order obtain the 
simple form shown. However, the examination the omitted 
terms made the appendix indicates that nearly all practical 
cases their omission justified. 


Two different types heat effects appear these equations. 
The first the change the steady-state liquid and vapor flow 
rates and (L; upon passage through the plate. 
This the familiar steady-state behavior usually resulting from 
differing heats vaporization among the plate components. 


The second effect, denoted the function jw), wholly 


transient phenomenon. Physically, Q(jw) associated with 
transient boiling and condensing action tray. When the 
concentration one component constant pressure boiling 
mixture fluctuates, the temperature the liquid the tray 
must also fluctuate the vapor pressure that pure component 
differs from the vapor pressure the mixture the mean 
plate temperature. Since the mixture boiling, the fluctuating 
temperature causes flashing and condensing the system which, 
turn, fluctuates the flow rates liquid and vapor. The heat 
source sink the sensible heat the mass the liquid 
the plate and the plate metal. 


These simulation results have implications the field 
distillation column control. common practice assume 
that liquid and vapor flow rates the column vary only when 
reflux and reboil conditions are varied. The present results 
indicate composition fluctuations can also induce considerable 
variation flow rates and this possibility should not ignored. 
Frequency response tests small laboratory size plate which 
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well mixed can used conjunction with Equations (24) 


and (25) determine the value and hence the import- 
ance this effect. 

The response the simulated extractive distillation tray 
was used with Equation (26) predict the plate efficiency, 
and this value turn was used obtain the equilibrium curve 
means Equation (20). The prediction proved excel- 
lent. 

The energy effects test used here unusually severe 
one. 90% the liquid the plate was high boiling extractive 
agent and slight fluctuations plate temperature caused rela- 
tively large quantities hydrocarbons flash condense 
order supply heat remove heat from the large mass 
low volatility extractive agent. ordinary distillation 
system where all the components exhibit comparable volatil- 
ties, the effects would much smaller and many systems 
insignificant. 


Energy Effects Imperfectly Mixed Plate 


Unfortunately, the problem treating the energy effects 
plate which imperfectly mixed has not yet been resolved 
satisfactory manner. Until such solution obtained, the 
frequency response technique limited perfectly mixed 
plates, small laboratory sized plates, systems 
which the heat effects are small, gas absorption plates and 
distillation systems which the volatilities and latent heats 
all components are comparable magnitude. Despite these 
limitations, the method its present form would appear 
adjunct the standard procedures. Further analytical and 
experimental work indicated. time hoped that 
frequency reponse procedure general validity can developed 
and tested. 
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Nomenclature 


amplitude ratio the effluent liquid concentration 
sinusoid the feed liquid concentration sinusoid, 
concentrations expressed mol fractions 


amplitude ratio the effluent liquid concentration 
sinusoid the feed vapor concentration sinusoid, 
concentrations expressed mol fractions 

constant 

molal heat capacity 

heat capacity plate which thermally interacts with 
fluids the plate 

CL = Crk + Cz2(1 

AC, = Cn — Cr 

ACy = Cy — 

point efficiency 

molai enthalpy 

molal flow rate liquid 

liquid level 

slope linear approximation the vapor-liquid 
equilibrium line 

the mixing parameter Kirschbaum pool concept: 
the total number equally sized, perfectly mixed 
stages which represents the liquid mixing the plate 

universal gas constant 

Laplace transform variable 

absolute temperature liquid 

temperature 

molal flow rate vapor 

liquid mol fraction 

vapor mol fraction 

vapor mol fraction equilibrium with liquid 

length liquid flow path contact with vapor 
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Greek Letters 


phase angle 

constant. Each value (1,2,...) denotes 
different constant. 

activity coefficient 

time 

molal heat vaporization 

= My — v2 

total pressure 

sinusoidal angular frequency 

Subscripts 

index liquid mixing stages, counting from the inlet 

liquid 

plate metal 

partial molal quantity 

vapor 

Superscripts 

time mean value 

time and area mean value 

the time-dependent component variable 

reference condition 

constant 
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APPENDIX 


Generalized Linear Analysis Binary, 
Perfectly Mixed Plate 


The objective the analysis specify the dynamic cause- 
and-effect relationships amongst the variables the plate. 
simultaneous set algebraic, transcendental, and ordinary dif- 
ferential equations written which approximates the 
behavior and implies the cause-and-effect relationships. When 
the set written sufficient detail include most the pre- 
dictable interactions, its complexity such that exact analytical 
solutions are generally impossible. Approximate analytical solu- 
tions may obtained, however, the equations are linearized 
time about steady-state operating condition. Since the linear- 
ization valid only the system makes small excursions from 
the steady-state, the approximate solutions are also valid under 
the same restriction. 


The cause-and-effect relationships amongst the system’s vari- 
ables are sought terms the transfer functions complex 
gains. The method involves the writing the approximate 
linearized set equations from the original set, the Laplace 
transformation the linear equations, and finally the rearrange- 
ment the transformed equations the form transfer 
functions between the desired pairs input-output variables. 
found, however, that when system involves fair number 
interacting variables, the practical problem rearranging 
manipulating the transformed equations into suitable transfer 
functions not easily accomplished not attacked systematically. 
order facilitate these manipulations, procedure devised 
Professor employed. The technique, called 
flow reduces set rules for the manipulation 
transformed, linear equations that desired transfer functions 
may extracted. the following analysis, the differential 
equations are written, linearized, and transformed. Then the 
flow graph technique explained and finally applied for the 
solution the three desired transfer functions. 

Rather than proceed writing the entire set equations, 
linearizing the entire set, and then transforming the entire set, 
these three operations are performed each individual equation 
turn. the original equations, each variable function 
time unmarked with either bar prime. When the 
equations are linearized, each variable replaced the sum 
its mean value denoted bar and its time-dependent compo- 
nent denoted prime. When nonlinearity occurs which 
involves the product dependent vaiables and those variables 
are replaced the sum their mean and time-dependent 
components, the expansion the product includes products 
time-dependent components. The linearization effected 
neglecting such terms since for small excursions, these are small 
comparison the other terms the equation. Each equation 
then translated zero mean subtracting from the 
steady-state form the original equation. The result subjected 
the Laplace transformation which yields equation relating 
the transformed, time-dependent components the original 
variables. These transformed variables are functions the in- 
dependent complex variable and are denoted simply the 
functional symbol For the transformation differential 
equations, the steady-state operation the plate about which 
the linearization effected represents convenient set initial 
conditions. Physically, this means that time the 
sinusoidal disturbance initially introduced the plate which 
has been running quiescently. Mathematically, specifies 
zero initial value for the time-dependent components each 
variable. 


Overall mass balance: the vapor holdup insignificant 
relative the liquid holdup, then 


Subtraction the steady-state form leaves 


which terms time varying components. With the initial 
value H’(0+) Laplace transformation yields 


Component balance: Writing mass balance for either 
the two species present, and assuming perfect mixing the 
liquid, results 


Substituting for each variable the sum its mean and deviating 
components 


dH’ 


Neglecting products primed quantities, the time-dependent 
results 


H— = x; L;+%,L; +9; +9; Vi —x'L—FV'—y’ 


Transformation with respect quiescent initial conditions yields 


Enthalpy balance: addition perfect mixing, the 
enthalpy balance assumes: (1) adiabatic operation, (2) the plate 
metal temperature equal that the liquid the plate, and 
(3) negligible variation vapor holdup enthalpy comparison 
the total enthalpy flux. 


Considering constant and substituting variable components 
yields 


The linearization 


The transformed relation can written 


Molal enthalpy expressions: For each the two compo- 
nents, zero enthalpy assigned the pure liquids base 
temperature 


Another pair analogous equations representing and hy; 
can written which all the variables are subscripted with 
the heat capacities are assumed independent temperature, 
the linearized time-dependent components the molal enthalpy 
expressions are given 
hy! [Cr + Cz2(1 x)] + (Ci Cre) x’ (7 t°) + 

Rh’ + (1 x) h' m2 + (Ami ie 


and two similar expressions for and 


Transformation and definition new terms (see Nomenclature) 
yields 


(s) and hy; (s) are written the same form Equations 
(A-4) and (A-5) but with subscripts all relevant symbols. 


Vapor-liquid equilibrium: linear vapor-liquid equilib- 
rium written: 


Subtracting the mean values and transforming the result yields 


Plate efficiency: The expression relating and 
involves constant plate efficiency value 


Transforming the time-dependent component relation 


Partial molal heats mixing and activity coefficients: 
the solutions are assumed then there excess 


The Canadian Journal Chemical Engineering, February, 1961 


pa rtia! 
mixing 


The 
and te 


where 
heat 


The 


The 
explic 
into 


smalle 
the 


Trans 


total 


The 


sume 


Tran 


The 


partial molal entropy for each component and the heats 
mixing are given 


The simplest correlation activity coefficient with composition 
and temperature Margules-type equation 


where and are constants. Eliminating from the 
heat mixing equations, and then linearizing 


The transformed heats mixing are 


The transformed activity coefficients must also expressed 
explicitly. Substituting mean and time-dependent components 
into the activity coefficient expressions, there results 


Since pseudo-linear behavior requires each much 
smaller than its the logarithms y’) are approximately 
Making these substitutions and proceeding with 
the linearization 


+ t' in = — 2a, Rx’ 


Y'2 + t' ln = 2a2 Rx’ 


Transforming 
1 


Pressure-temperature-composition relationship: The 
total pressure the sum the two partial pressures. 


= yiPix + 2P2 (1 — x) 


The transformed relation 


Vapor pressure: The vapor pressure each species as- 
sumed satisfy the Clausius-Clapeyron approximation. 


In Py = as — im Ps 
Proceeding with the linearization 
= Ay; 
F; P; = 
n(P; + i’) R(T + T’) 
Ay; 
>; R T? 
Transforming 


Pressure drop: The pressure drop suffered the vapor 
passes through the plate assumed depend both upon the 
liquid level the plate and upon the vapor flow rate. 


Tw; =a; +a; V;? 
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The linearization 
The transform the linearization 
2a; 


Liquid level: The liquid level assumed depend upon 
the liquid flow rate some positive power 
This equation can approximated linearly 


which transforms into 


holdup: Liquid holdup related linearly liquid 
eve 
H= + Qiol 
H' = Gio I’ 


One additional interaction can predicted but has not been 
included. That the relationship between the pressure under the 
plate and the incoming vapor flow rate V;. Clearly, variations 
liquid level must introduce variations pressure drop across 
the plate and pressure drop fluctuations must reflected 
either the vapor throughput rate the under-plate pressure 
both. However, the form this effect depends almost 
entirely upon the characteristics the vapor supply equipment. 
all likelihood the interaction will not sufficiently strong 
influence the dynamic behavior the plate any great extent 
but should kept mind when plans are drawn for 
experimental vapor supply system. 

The transformed equations governing the linearized behavior 
the distillation plate comprise Equations (A-1) through (A-14). 
this point, convenient eliminate several the output 
variables and thereby reduce the number equations. (s) 
may eliminated from Equations (A-6) and (A-7) that 


The heats mixing are eliminated substitution Equations 
(A-8) into (A-4) 
=> Cr t(s) (AC,At + Ahm) x(s) 


Vapor pressures, Equation (A-11) and activity coefficients, Equa- 
tion (A-9) may substituted Equation (A-10) 


The liquid level is.substituted both Equations (A-12) 
and (A-14) 


A(s) = AsQio Lis) (A-19) 


The remaining transformed variables are now categorized 
either input output. The input variables are, sense, 
independent since they not depend upon the dynamic behavior 
the plate. The output response variables, the other 
hand, must completely specified the set equations. 

The input variables are V;, x;, and The 
output variables are hz, hy, and These are 
specified Equations (A-19), (A-3), (A-1), (A-2), (A-15), (A-16), 
(A-5), (A-17) and (A-18) respectively. 


this analysis, x;, and hy; are the only input 
variables which disturbances are introduced the system. 
The others are held constant and hence their time-deviating 
components and their transforms are zero. Consequently, 

the flow graph method for manipulating the equations, 
the interrelationships amongst the variables are represented 
network directed branches which connect nodes. The 
nodes represent the variables; each branch represents the cause- 
and-effect relationship gain with respect the nodes 
connects. One starts the flow graph establishing node for 
each variable the set equations, both input and output 
variety. The system defined nine equations, one for each 
output variable. example, Equation (A-1) specifies the 
variable H(s). says that each the variables L(s) 
and V(s) affects H(s). ignored.) Consequently, three 
directed branches are drawn the flow graph from L(s) 
and V(s) H(s). The gains associated with these branches are 
given the equation. For example, the gain the L(s) H(s) 
branch 1/s) since this when multiplied expresses 
the L(s) H(s) relationship the equation. the same manner, 
all the relationships are entered branches the flow graph. 
This shown for the entire set equations Flow Graph 
Figure The branch gains are listed the format: 


1/s which specifies gain 1/s for the 
branch from node L(s) node H(s) Flow Graph 


From Equation (A-1) which specifies H(s): 


1/s 


From Equation (A-2) which specifies x(s): 


From Equation (A-3) which specifies V(s): 

H(s)—> V(s) = h,/hy 


V/hy 


From Equation (A-5) which specifies hy(s) and 


From Equation (A-15) which specifies 


From Equation (A-16) which specifies and 


From Equation (A-17) which specifies 


From Equation (A-18) which specifies 


From Equation (A-19) which specifies L(s): 


Flow Graph which all the above branches are drawn, 
topologically represents the set transformed equations. The 
first step the reduction such graph the formation 
new flow graph containing fewer nodes than the first. The 
nodes which are retained, called residual nodes, are the sources, 
the sinks, and minimum set nodes whose members participate 
all the closed loops the first graph. source node 
into which arrows enter, and sink one from which 
arrows emanate. closed loop series branches connecting 
two more nodes which may travelled around endlessly, 
e.g., branches connecting L(s) V(s) H(s) 
L(s). examination Flow Graph shows that sources 
are and y,(s), that L*(s) the only sink, and that 
least either H(s) x(s) part every closed loop. result, 
all the nodes except H(s) and x(s) can elim- 
inated. the elimination node corresponds 
that the retention the minimum five nodes just listed means 
the elimination But constitutes output variable 
which must appear one the final transfer functions. 
order save L(s), new variable L*(s) introduced which 
related L(s) L*(s) L(s). L*(s) node entered 
Flow Graph and connected (Z(s) branch whose gain 
unity. 


Clearly, sink and residual node. 

The residual nodes H(s), x(s) and are 
entered Flow Graph branch appears Flow Graph 
connecting two residual nodes if, and only if, Flow Graph has 
one more residual paths connecting the same two nodes. 
residual path defined connected series branches running 
from residual node itself another residual node but 
not passing through any other residual node. 


order determine the gains the branches Flow 
Graph only two transformations are needed: 


4) g3 


cascade transformation (a), the resultant gain the 
product the component gains. the multipath transformation 
(b), the resultant gain the sum the component gains. 


With these transformations, the branches Flow Graph 
can calculated from Flow Graph branches. The following 
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first tc 
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lis: shows the Flow Graph branch gains terms Flow Graph 
gains: The functional notation suppressed for convenience. 


(A-21) 


These transfer functions represent the overall cause-and-effect 
relationships between the indicated input and output variables. 
Included part each transfer function that simplified 
cause-and-effect relationship which forms the basis the dynamic 
analysis method. The remainder each overall transfer function 
represents dynamic behavior which interferes with the analysis 
method and, therefore, must small for the analysis 
successful. For example, Equation (A-20) describes the overall 
effect the liquid effluent composition when the incoming liquid 
composition disturbed. its present form, Equation (A-20) 
composed three groups interactions, the two terms 
the numerator and the second term the denominator. only 
the first term the numerator expanded terms distillation 
parameters, one sees that the desirable part the overall transfer 
function included this group. 


Flow Graph has sufficiently simple form that the desired 
transfer functions may obtained inspection. necessary 
first mention one more transformation: 


1 — ge 1 — go 


node with self-loop gain may replaced cascade 
branch gain With this self-loop transformation, the 
desired transfer functions terms Flow Graph gains are: 
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The expression enclosed dashed lines that part the overall 
transfer which the basis the dynamic analysis method (see 
Equation (9)). Clearly, then, all the remaining parts Equatio 
(A-23) plus the two groups terms 


Equation (A-20) represent undesirable interactions. Therefore, 
the dynamic analysis method unharmed extraneous 
effects, 


must good approximation Equation (A-20). This will 
the case 


and 


Equation (A-21) analyzed this same manner, 
seen that for the desirable part 


good approximation the overall transfer function, the 
inequalities (A-25) and (A-27) and also 


and 


must satisfied. 
Some comments can made the nature these inequal- 
ities: The expressions 


and 


the inequalities (A-26), (A-27) and (A-29), respectively, repre- 
sent the dynamic behavior which can pictured crudely 
follows: The fluctuation composition y;) causes 
fluctuation the number mols holdup the plate; then 
the reaction the holdup fluctuation fluctuation the 
effluent liquid composition. intuitive reasoning, neither 
these two effects should strong. The molal holdup 
distillation plate depends primarily flow rates but also 
the relative enthalpies the feed and effluent streams. The 
stream enthalpies, course, depend upon composition, but when 
the component species are thermodynamically similar, this de- 
pendency weak. highly unlikely that reasonable composi- 
tion fluctuations propagate strongly the holdup. 
reasoning correct, then the inequalities (A-26), (A-27) and 
(A-29) should easily satisfied. 

inequality (A-24) written terms distillation para- 
meters, takes the form 


This condition should also satisfied since the right-hand 
side made small the large denominator the same 
reasoning, inequality (A-28) should satisfied cases. 
can written 


The remaining inequality (A-25). Since has been reasoned 
that all the other inequalities should easily satisfied, the 
principal liquid and vapor fluctuation effects must felt here. 
this is, indeed, the case, then fortunate situation occurs: 
Equations (A-20) and (A-21) reduce 


x(s) 


and 


where constitutes all the terms representing the liquid and 
vapor effects and the right hand side the in- 
equality (A-25) 

The ratio equations 


A-32 
x(s)/x,(s) L; (AR);;: 


Equations (A-20), (A-31) and (A-32) with replaced are 
equations (24), (25) and (26) the body this article. 
formally the terms the right-hand side the inequality (A-25). 
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The Horizontal Pipeline Flow Equal 
Density 


CHARLES?, 


and 


The horizontal flow equal density oil-water 
mixtures was investigated l-inch diameter 
laboratory pipeline. Oils viscosities 6.29, 16.8 and 
65.0 centipoise were used the experiments. Flow 
patterns, holdup ratios and pressure gradients were 
investigated for range superficial oil velocity 
from 0.05 3.0 and range superficial 
water velocity from 0.1 3.5 ft./sec., with input 
oil-water ratios ranging from 0.1 10.0. 


Similar series flow patterns were observed for 
each oil and were found largely independent 
the oil viscosity. high oil-water ratios oil formed 
the continuous phase and water-drops-in-oil regime 
was noted. the oil-water ratio was decreased the 
flow patterns concentric oil-in-water, oil-slugs-in-water, 
oil-bubbles-in-water and oil-drops-in-water which 
water was the continuous phase were observed. The 
most viscous oil did, however, exhibit anomalous 
behavior low superficial water velocities and this 
attributed different interfacial properties. 


The pressure drops measured indicated that for 
given oil flow rate the pressure gradient was 
reduced minimum the addition water 
provided that the oil was not turbulent flow. 
Maximum pressure gradient reduction factors varied 
from 1.7 for the 6.29 centipoise viscosity oil 
for the 65.0 centipoise viscosity oil input oil-water 
ratios 4.5 and 1.0 respectively. 


contact two immiscible liquids encountered widely 
the chemical and petroleum industries. The design 
tubular reactors, liquid-liquid extraction equipment and pipelines 
for liquid transportation may depend knowledge the 
flow characteristics the simultaneous movement immiscible 
liquids. Although there considerable amount published 
information the flow liquid-gas and liquid-solid mixtures 
there very little the flow liquid-liquid mixtures. 


knowledge the factors affecting the transport crude 
oil and water pipelines important the petroleum industry. 
Water often produced large quantities with crude oil and 
the characteristics this two-phase flow are interest both 
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the well itself, where the flow vertical the production 
tubing, and horizontal pipelines transporting the crude oil 
field treating facilities. The major interest however has been 
concerned with the beneficial effect, way decrease 
pressure drop, which results from the introduction water 
controlled amounts into pipelines carrying heavy viscous crude 
oil. The introduction water can reduce the gradient 
along the and the power requirement necessary 
pump given quantity oil. 

The flow pattern which appears most attractive from the 
viewpoint pressure gradient reduction the addition 
water viscous crude oil that the concentric flow oil 
and water which the less viscous water forms uniform 
annulus the region high shear rate next the pipe wall. 
Russell and made theoretical study the laminar 
concentric flow two liquids equal densities but different 
viscosities smooth circular pipe and also the laminar 
stratified flow between wide parallel plates. Concentric flow 
was also analyzed Chernikin® for the case when the water 
turbulent flow and the oil laminar flow. Patents 
Clark and and Chilton and deal with the 
commercial application water reduce oil pipeline pressure 
gradients. and were both concerned with the 
relative permeability oil and water porous media and made 
theoretical analyses the flow idealised capillary systems. 
measured pressure drops and heat transfer the 
pipeline flow oil-water emulsions. 


general, two immiscible liquids different densities tend 
stratify when flowing together pipeline. The greater the 
density difference the more nearly complete the stratification. 
the case oil and water, the oil usually has density less 
than that the water and forms upper layer flowing above 
the water. Stratified flow has been investigated experimentally 
Russell, Hodgson and with light oil horizontal 
pipeline one-inch diameter and with heavy 
crude oil and pipelines. Charles and 
numerically analysed stratified flow circular pipe with the 
aid digital computer. 

The relationship between the pressure gradient and flow rate 
for single liquid usually shown plot friction factor 
against Rey nolds number. These quantities incorporate the 
pressure gradient, the liquid velocity, the pipe diameter, and 
the density and viscosity the liquid. The friction factor 
and the Reynolds number cannot used the normal way 
for the simultaneous flow two liquids. Also, the interfacial 
tension the two liquids, which the ease emulsifica- 
tion, and the contact angle the liquid-liquid interface the 
tubing material, which determines the liquid that preferentially 
wets the pipe, may important variables, particularly low 
liquid flow rates when the surface forces are longer negligible 
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Figure 1—Schematic diagram the pipeline and accessory equipment used the study the simultaneous flow oil 
and water equal densities. 


compared with the total kinetic energy. the present investi- 
gation the densities the liquids were equal and there 
difficulty representing the density the two expressions. 
However, difficulty arises representing the velocity and 
viscosity terms the composite flowing stream. Russell 
investigation the flow oil and water with 
different densities presented their pressure drop data terms 
friction factor calculated the basis the superficial 
velocity the water and the density the water. This tech- 
nique has correlating advantages but does not permit ready 
visualization the effect the introduction water the 
pressure gradient flowing oil stream. The more direct 
relationship pressure gradient the ratio the phases 
constant superficial velocity one them useful demon- 
strating the practical effect. This method used reporting 
the results the present investigation. 


While the word “superficial” used modify the oil and 
water velocities will omitted throughout the remainder 
this discussion, all liquid velocities referred will superficial 
the sense that they are calculated from the flow rate and total 
cross sectional area the pipe. 


Other quantities which are convenient employ the 
treatment two-phase flow are the input ratio, the situ 
ratio and the holdup ratio. The input oil-water ratio, 
defined the ratio the oil flow rate the water flow rate. 
The situ ratio the ratio the volume the pipe occupied 
oil the volume occupied water. general the input 
and situ ratios differ because the liquid contact with the 
pipe wall will tend accumulate the pipe, and the ratio the 
input oil-water ratio the situ oil-water ratio defined 
the holdup ratio. Thus, general, for given input ratio, 
the oil wets the pipe the holdup ratio will tend less than 
unity because the situ ratio greater than the input ratio and, 
the water wets the pipe the holdup ratio will tend greater 


than unity because the situ ratio less than the input ratio. 
will also convenient consider the input water-oil ratio, 
which the reciprocal the input oil-water ratio. 


Experimental Equipment and Procedure 
Test Oils 


Three Newtonian oils were used the investigation and 
all were obtained increasing the density three commercial 
oils (Marcol GX, Wyrol and Teresso 85, marketed 
Imperial Oil Ltd.). Sufficient carbon tetrachloride was added 
each oil give the oils the same density the water (0.998 
gm. per temperature 25°C.). Carbon tetrachloride, 
with density 1.59 gm. per cm.* was selected the basis 
high density and relatively low cost the most suitable liquid 
increase the density the oil phase. Precautions were taken 
confine the toxic carbon tetrachloride vapor protect the 
operating personnel. The properties the test oils are given 
Table 


Equipment 

flow diagram the pipeline and accessory equipment 
shown Figure test oil and water supplies were 
contained storage tanks and the oil was maintained 
uniform condition means double propellor-type mixer. 
The oil temperature was adjusted recycling the oil through 
heater, until the temperature was indicated 25°C. 
bimetallic thermometer situated the recycle line. The 
water temperature was controlled adjusting the relative 
amounts hot and cold water flowing into the storage tank. 
The density the oil was checked using hydrometer 
sample taken from the storage tank and was usually found 
slightly low (for example, 0.995 gm. per compared 
with the required value 0.998 gm. per cm.*) because the 
loss carbon tetrachloride through evaporation. Sufficient 


TABLE 


PROPERITES 


Approximate 

3 
cp. 25°C. 25°C. content, volume per cent oil used base 
6.29 0.998 clear 20.6 Marcol 
16.8 0.998 clear 18.7 Wyrol 
65.0 0.998 dark green-brown 16.7 Teresso 
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tetrachloride was added give the oil supply the re- 
density. 

flow oil and water the pipeline was established 
screw-type positive displacement pumps which delivered the 
and water separately calibrated rotameters after which 
the oil and water flows were combined concentrically the 
beginning ten-foot long section means nozzle 
which the oil was introduced inside water annulus. 


The pipeline calming section and twenty-four foot test 
section were supported horizontal position wooden 
trestle. cellulose acetate-butyrate pipe was smooth and 
transparent with internal diameter 1.04 in. The pipeline 
was mounted brass slides and with flexible coupling 
between the pipeline and the separation unit that was free 
expand under the influence ambient temperature. The 
test section was supplied with pressure taps for the measure- 
ment the pressure drop and with two quick-action valves 
joined taut steel cable for the measurement the situ 
contents the pipeline. 


The oil was separated from the water for return storage 
establishing density differential between the immiscible 
liquids the addition saturated brine the oil and water 
immediately after left the test section. The brine 
increased the density the aqueous phase which then formed 
the lower layer the separator. upper oil layer was 
pumped storage through salt tower which removed the 
last traces water present small drops. 


Observation flow pattern 


The oil and water flow rates were adjusted the desired 
values and time was permitted for steady conditions 
established usually one five minutes. general, the slower 
the flow rates, the longer was the time required 
steady conditions. When visual observation indicated steady 
conditions the flow pattern was noted. Photographs were 
taken representative flow patterns. 


Pressure drop measurement 


The pressure drop measuring system was similar that 
used Russell and consisted two pressure stations, 
situated 24.0 feet apart the test section, with copper connect- 
ing lines manometers. The fluid the connecting lines was 
air, the pressure which was adjusted pressure reduction 
valves connected per gauge air supply. The oil 
and water mixture was admitted into the splash pot each 
pressure station and the liquid levels adjusted the same 
elevation manipulating the air pressure reduction valves 
until the connecting line pressure balanced the pipeline pressure. 
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Figure Comparison theoretical and experimental 
friction factors for the flow water and the 6.29, 16.8 
and 65.0 centipoise viscosity oils. 
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The pressure differential between the pressure taps was indi- 
cated water mercury manometer. 


Measurement situ contents 


After closing the pressure taps the quick action valves 
the ends the test section were closed simultaneously and the 
oil and water pumps stopped. which fitted 
the internal surface the pipeline, was introduced into the 
pipeline upstream the test section, and, when forced through 
the test section air pressure, ejected the trapped oil and ater 
into volume measuring equipment. The passage the 
provided fast, efficient method ejecting the oil and water. 


Interfacial tension measurement 


The oil-water interfacial tensions were determined using 
simple capillary which involved the measurement 
the equilibrium heights oil and water columns having joint 
interface inside glass capillary known diameter. 


Experimental Results and Discussion 


The operation the pipeline, including the rotameter 
calibrations and the manometer system for measuring pressure 
drops, was initially checked recording pressure “drops for 
the flow water and the three test oils separately. Fanning 
friction factors were plotted against Reynolds numbers 
logarithmic paper. The good agreement between the experi- 
mental data and the accepted curves for the laminar and turbulent 
regions evident Figure 

was desired observe the flow the oils with water 
over wide ranges oil and water flow rates possible and 
obtain pressure drop and holdup measurements for the various 
flow patterns observed. Flow patterns were observed over 
range oil velocity from 0.05 3.0 ft./sec. and range 
water velocity from 0.1 3.5 ft./sec. Pressure drop and holdup 
measurements were made about five oil velocities for each 
five water velocities. 


Flow patterns 


series flow patterns was observed for decrease 
the oil flow rate constant water flow rate. The flow pattern 
changed from water oil, through concentric, 
oil-slugs-in-water and oil-bubbles-in-water, oil-drops-in-water. 
Drawings were prepared from photographs taken for the 16.8 cp. 
viscosity oil, and are representative for the three oils except for 
the flow patterns the most viscous oil low water flow rates. 
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Figure Drawings prepared from photographs the 

16.8 viscosity oil flowing the presence water and 

showing the variation flow pattern with oil velocity for 
low fixed water velocity 0.10 ft./sec. 
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Figure Drawings prepared from photographs the 

16.8 centipoise viscosity oil flowing the presence 

water and showing the variation flow pattern with oil 
velocity for fixed water velocity 0.682 ft./sec. 
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Figure Drawings prepared from photographs the 

16.8 centipoise viscosity oil flowing the presence 

water and showing the variation flow pattern with oil 
velocity for high fixed water velocity 2.04 ft./sec. 
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Figure 6—Flow regimes for oils viscosity 6.29 and 16.8 
centipoise and density 0.998 flowing 
presence water with the same density. 


The Canadian Journal Chemical Engineering, February, 1961 


The effect the flow pattern reducing the oil velocity 
the relatively low water velocity 0.1 ft. per sec. shown 
Figure high oil velocities the water-drops-in-oil flow 
pattern was observed the example for oil velocity 
1.36 ft./sec. Reduction the oil velocity about 1.2 
caused both oil and water present continuous 
Concentric flow, which the oil formed continuous core 
along the axis the pipe and the water uniform layer the 
pipe wall, shown for oil velocity 0.454 ft./sec. 
further reduction oil velocity water velocity 0.1 ft./sec. 
resulted the oil phase becoming discontinuous. The oil 
flowed slugs, and lower oil velocities, bubbles. The 
oil-slugs-in-water flow pattern shown for oil velocity 
0.200 ft./sec. and the oil-bubbles-in-water flow pattern shown 
for oil velocity 0.050 ft./sec. 


The effect the flow pattern reducing the oil velocity 
higher water velocity similar, shown Figure for 
water velocity 0.682 ft./sec. Again high oil velocities 
the water was dispersed small drops and example this 
flow pattern given for oil velocity 1.95 ft./sec. 
reduction the oil velocity again resulted concentric oil-in- 
water flow and example shown for oil velocity 
1.11 ft./sec. The flow was disturbed compared with the con- 
centric flow water velocity 0.10 ft./sec.; the oil-water 
interface was wavy and some water was entrained the oil. 
further reduction oil velocity again resulted oil-slugs-in- 
water flow and subsequently flow. 
example slug flow shown for oil velocity 0.682 ft./sec. 
and example bubble flow for oil velocity 0.200 


The effect the flow pattern reducing the oil velocity 
even higher water velocity 2.04 ft./sec. shown 
Figure high oil velocities the water-drops-in-oil pattern 
was again encountered. the lower oil velocity 1.82 ft./sec. 
the considerably disturbed form concentric flow was encount- 
ered and further reduction oil velocity gave the oil-slugs-in- 
water, the oil-bubbles-in-water and the oil-drops-in-water flow 
patterns. this relatively high water velocity was difficult 
differentiate slug, bubble and drop flow patterns because, 
for given oil velocity, the oil was often present slugs, 
bubbles and drops various sizes shown for oil velocities 
0.682 and 0.200 ft./sec. The drawings indicate that for the 
concentric oil-in-water, oil-slugs-in-water, oil-bubbles-in-water 
and oil-drops-in-water flow patterns the pipe was always ina 
water-wet condition. 


The flow pattern regimes are summarized for the 6.29 and 
16.8 cp. viscosity oils Figure where the oil and water 
velocities are represented logarithmic coordinates and single 
lines are drawn between the various flow regimes. The change 
from one regime another was not well defined and the lines 
represent gradual transitions. The oil and water velocities for 
which pressure drop and holdup measurements were made, are 
plotted along with the corresponding experiment number 
Figure are also the velocities which are represented 
the flow pattern examples Figures and The flow 
regimes for these two oils are almost identical. The only 
difference occurs the position the boundary between the 
water-drops-in-oil and the concentric oil-in-water flow regimes. 
The boundary somewhat higher oil velocity for the 
16.8 cp. oil than for the 6.29 cp. oil. viscosity differential 
would appear not account for this difference since this bound- 
ary for the most viscous oil not displaced further but 
approximately the same position the boundary for the 6.29 cp. 
oil shown Figure 


The flow regimes observed for the 65.0 cp. oil are sum- 
marized Figure and show general agreement with those 
for the other less viscous oils except low water velccities 
where succession different flow patterns was obtained with 
the oil the continuous phase. The line represents the 
conditions for which the flow pattern was unstable and repre- 
sents inversion between flow patterns which the oil was 
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Figure 7—Flow regimes for oil viscosity 65.0 centipoise 
and density 0.998 flowing the presence 
water with the same density. 


the continuous phase and flow patterns which water was the 
continuous phase. For given oil velocity, water velocity 
less than that indicated the inversion line gave flow pattern 
with oil the continuous phase, while for water velocity 
greater than that indicated the inversion line, flow patterns 
with water the continuous phase were obtained. The bound- 
aries between the regimes which water was the continuous 
phase were almost identical with the corresponding boundaries 
for the less viscous oils. The regime boundaries are represented 
broken lines low oil and water velocities because the flow 
patterns were less well defined under these conditions. 


The flow patterns observed this study may related 
the flow patterns reported Russell examining the 
flow oil and water different densities: 


(a) bubble flow, which oil drops and bubbles moved 
along the top the pipe. 

(b) stratified flow, which the oil and water flowed 
two layers with quiet interface. 


(c) mixed flow, which occurred high flow rates and was 
completely mixed type flow. 


The bubble flow Russell appears correspond the 
water flow patterns observed the present study. The stratified 
flow corresponds the concentric oil-in-water flow pattern and 
the mixed flow the water-drops-in-oil flow pattern. 


circulation oil within the bubbles and slugs was noted 
for the oil-slugs-in-water and oil-bubbles-in-water flow patterns. 
This movement indicated Figure for bubble flow. 
Approximate velocities relative the average flow velocity 
are represented the vector arrows. The direction rotation 
the bubble probably dependent the direction small 
density differential between the oil and the water. When the 
oil slightly more dense than the water the bubbles would 
tend drag more the lower surface the pipe and rotate 
the direction shown. circulatory movement was also 
noticed the water phase between the bubbles slugs oil. 
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Figure Liquid movement the oil-bubbles-in-water 
flow pattern. 


Interfacial properties and their influence flow pattern 


has been observed this discussion that the most viscous 
oil showed behavior different from the less viscous oils low 
water velocities. Since under these conditions surface forces 
may longer negligible compared with the total kinetic 
energy the flow the dissimilar behavior the most viscous 
oil suggested that interfacial properties might important 
the determination flow pattern. The interfacial oil-water 
tensions (accurate within about five per cent) were measured 
and are reported Table The similar values for the 6.29 
and the 16.8 cp. viscosity oils and the dissimilar value for the 
65.0 cp. oil are not surprising, since the 6.29 cp. and the 16.8 cp. 
viscosity oils were both prepared from white mineral oils while 
the 65.0 cp. oil was prepared from lubricating oil containing 
anti-corrosion compounds. The interfacial contact angles with 
the pipe surface were observed qualitatively: the pipeline flow 
was stopped with water the continuous phase and the most 
viscous oil was visually observed make contact with, and 
spread readily on, the pipe surface, whereas the less viscous 
oils maintained definite angle contact with the pipe surface 
the presence water. The influence the different inter- 
facial tensions the flow patterns not clear since would 
expected that the most viscous oil having relatively low 
oil-water interfacial tension would exhibit greater tendency 
form emulsions. However the series flow patterns 
which the most viscous oil formed the continuous phase 
probably accounted for the preferential wetting the pipe 
surface this oil. 


Holdup ratio 

When two fluids flow together pipeline the situ 
volumetric ratio is, general, different from the input volu- 
metric ratio since the liquid contact with the pipe wall tends 
held back. For highly turbulent flow, however, with 
both fluids thoroughly mixed together, the situ volumetric 
ratio would approximately equal the input volumetric 
ratio. The holdup ratio has been defined the ratio the input 
oil-water ratio the situ oil-water ratio. The holdup ratio, 
then, does general differ from unity, and will tend 


TABLE 
INTERFACIAL TENSIONS 


Oil viscosity Interfacial oil-water tension 


cp. 25°C. dynes/cm. 25°C. 
6.29 
16.8 
65.0 


= 8— in water x K x Yy inte 
| in water 44 
a (| Yx xxx Vj 


TABLE (Partial) 


EXPERIMENTAL DATA FOR THE THE 16.8 AND WATER CIRCULAR PIPE 


Liquid Temperatures: 25°C. 
Pipe diameter: 1.04 in. 


Liquid densities: 0.998 gm. per 


Oil Viscosity: 16.8 cp. 


Water viscosity: 0.894 cp. 


Distance between pressure tappings: 24.0 ft. 


Volume pipeline for situ measurements: 4720 


Superficial Superficial Input oil- Pressure drop situ 
water velocity oil velocity water ratio over test section ratio 
ft./sec. ft./sec. ft.3/ft.3 in. water ft.3/ft.3 
0.100 1.00 10.0 2.50 6.68 1.50 
0.100 0.500 5.00 1.0-1.2 3.07 1.66 
0.100 0.200 2.00 0.45 1.56 1.28 
0.100 0.100 1.00 0.30 0.958 1.05 
0.100 0.050 0.500 0.30 0.428 
0.300 3.00 10.0 17.6 0.826 
0.300 1.50 5.00 3.45 3.26 
0.300 0.600 2.00 0.70 1.55 1.29 
0.300 0.300 1.00 0.35 0.775 1.29 
0.300 0.150 0.500 0.35 0.426 
0.800 3.98 4.98 21.0 3.97 1.26 
0.800 1.60 2.00 7.20 1.48 1.35 
0.800 0.80 1.00 2.55 0.840 1.19 
0.800 0.40 0.500 1.60 0.463 1.08 
0.800 0.159 0.199 1.10 0.178 
1.80 3.60 2.00 39.2 1.70 1.18 
1.80 1.80 1.00 14.2 0.796 1.26 
1.80 0.900 0.500 8.70 0.440 1.14 
1.80 0.358 0.199 6.20 0.174 
1.80 0.182 0.101 5.20 0.094 1.06 
3.50 3.50 1.00 40.0 0.953 1.07 
3.50 0.500 27.3 0.407 
3.50 0.693 0.198 19.6 0.189 1.06 
3.50 0.350 0.100 0.094 1.07 


| 
| 
| 
| 
| 


greater than unity when water the component contact 
with the pipe wall and less than unity when oil 
contact with the pipe wall. 

Holdup ratios were measured various oil-water input 
ratios and water velocities 0.1, 0.3, 0.8, 1.8 and 3.5 ft./sec. 
for each the three oils. These water velocities were chosen 


1.8 

3.5 

4 

ft/sec 
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INPUT RATIO, 
Figure Holdup ratio function the input oil- 


water ratio and water velocity for oil viscosity 6.29 
centipoise and density 0.998 


give suitable separation the data logarithmic scale. 
example the experimental data obtained given 
Table where the holdup ratios and pressure drops for the 
16.8 cp. viscosity oil are presented*. The experimental holdup 
ratios are plotted against the input oil-water ratio, with water 
velocity parameter, for the 6.29 cp. viscosity oil Figure 


2.0 
=2.0ft/sec 


INPUT RATIO, 


Figure 10—Holdup ratio function the input oil- 
water ratio and oil velocity for oil viscosity 6.29 centi- 
poise and density 0.998 
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Figure 11—Holdup ratio function the input oil- 
water ratio and oil velocity for oil viscosity 16.8 centi- 
poise and density 0.998 


(In Figures 13, and the solid data points represent 
experiments which the oil wetted the pipe wall and the open 
points represent experiments which the water wetted the 
pipe.) The holdup curves pass through maxima, the locations 
which coincide with the concentric oil-in-water and oil-slugs- 
in-water flow patterns. high flow velocities the holdup 
ratios tend become independent the flow rate either 
water. This more clearly indicated Figure 10, where 
the holdup ratio plotted against the input oil-water ratio with 
oil-velocity parameter. Figure was prepared from the 
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Figure 13—Experimental pressure gradients for the oil 
viscosity 6.29 centipoise function the input oil- 
water ratio and water velocity. 
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Figure ratio function the input oil- 
water ratio and oil velocity for oil viscosity 65.0 centi- 
poise and density 0.998 


smooth curves drawn for the constant water velocities 
Figure tendency for the holdup ratio approach 
unity low oil-water ratios well illustrated Figure 
for the range oil velocities investigated. 

Figures similar Figure were prepared from the plots 
the experimental data constant water velocity for the 16.8 
and 65.0 cp. viscosity oils. The dependence the holdup ratio 
oil velocity high oil-water ratios indicated for the 16.8 cp. 
oil Figure 11, with the tendency approach value 
unity low oil-water ratios. The holdup ratio values for the 
most viscous oil Figure reflected the different behaviour 
this oil high oil-water ratios. The holdup ratio values 
were less than unity for the flow regimes which the oil 
formed the continuous phase and showed little dependence 
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Figure 14—Pressure gradients for the oil viscosity 6.29 
centipoise function the input water-oil ratio and oil 
velocity. Points are calculated from the curves Figure 13. 
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Figure 15—-Experimental pressure gradients for the oil 
viscosity 16.8 centipoise function the input oil- 
water ratio and water velocity. 


the oil velocity. low oil-water ratios the holdup ratio values 
were very similar the values obtained for the less viscous 
oils the same oil-water ratios, and tended approach unity 
the oil-water ratio decreased. 


Pressure gradients 


Pressure drops were measured the same oil and water 
flow rates those for which holdup ratios were recorded 
various input oil-water ratios for water velocities 
0.1, 0.3, 0.8, 1.8 and 3.5 ft./sec. The pressure drops obtained 
for the 16.8 viscosity oil are given Table 3.* 


The pressure gradients obtained for the 6.29 cp. viscosity 
oil are presented Figure where the pressure gradient 
plotted against the input oil-water ratio logarithmic scales 
with water velocity the parameter. The points zero 
oil-water ratio represent conditions for which water flows alone 
and were calculated from the friction factor Reynolds number 
curves for the flow single liquid smooth circular pipe. 
The experimental data for water velocities 0.3, 0.8, 1.8 and 
3.5 conform family relatively smooth curves. 
The data for water velocity 0.1 ft./sec. are consistent 
with the data for higher water velocities only oil-water 
ratios above about for lower input ratios the pressure gradients 
observed appear three five times larger than those 
expected from extrapolation the family curves from 
higher water velocities. The pressure gradients water 
velocity 0.1 ft./sec. were some the lowest measured but 
there reason suspect the accuracy the 
possible that the apparently high pressure gradients may 
for the tendency, these low oil and water 


*Complete experimental data Table this paper has been deposited 
as Document No. 6438 with the ADI Auxiliary Publications Project, 
Photoduplication Service, Library of Congress, Washington 25, D.C. 
A copy may be obtained by citing the Document No. and by remitting 
$1.25 for photoprints, or $1.25 for 35 mm. microfilm. Advance payment 
is required. Make cheques or money orders payable to: Chief, Photo- 
duplication Servce, Library Congress. 
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Figure 16—Pressure gradients for the oil viscosity 16.8 
centipoise function the input water-oil ratio and oil 
velocity. Points are from the curves Figure 15. 


flow rates, for the oil phase drag the pipe wall result 
slight density differential between the oil and water phases. 

Since primary interest lies the pressure gradients for 
given oil flow rate the presence various amounts water, 
the results are better plotted, Figure 14, the pressure 
gradient versus the water-oil input ratio with the oil velocity 
the parameter. The curves were constructed from the smooth 
curves Figure and the points not represent direct 
experimental observations. The points intersection the 
curves with the ordinate were calculated for oil viscosity 
6.29 cp. from the accepted friction factor Reynolds number 
curve. 


apparent that for oil flow rates below ft./sec. re- 
duction pressure gradient obtained the addition water. 
The factors which the pressure gradients are reduced are 
largest the lowest oil flow rates. The maximum factor 
which the pressure gradient was reduced was approximately 1.7 
oil velocity 0.5 ft./sec. and input water-oil ratio 
0.22 (or water percentage 18% based the total liquid 


For oil velocities above 1.0 ft./sec. pressure gradient 
reduction was not observed and increase pressure gradient 
was indicated immediately upon the addition water. Maxima 
occurred the curves input water-oil ratio 0.1 and are 
probably accounted for transition from laminar turbulent 
flow the addition the water. This likely so, 
especially oil velocity 1.5 ft./sec. since, this case, 
the Reynolds number for the the oil nearly 2000. 
high water-oil ratios the curves show steady increases 
pressure gradient with increase water-oil ratio and ratio 
the presence the oil had little effect the pressure 
gradient. That is, the pressure drops recorded and the calculated 
pressure drops assuming the liquid all water (rather than 
parts water and part oil) agree within few per cent. 

The flow pattern loci from Figure are superimposed 
Figure and evident that, the oil flow rates for which 
the pressure gradient was reduced below that for the oil flowing 
alone, the pressure gradient was least the concentric oil-in- 
water and oil-slugs-in-water patterns. 
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Figure 17—Experimenal pressure gradients for the oil 
viscosity 65.0 centipoise function the input oil- 
water ratio and water velocity. 


The data for the 16.8 cp. oil are shown Figure 15, where 
again the pressure gradient plotted against input oil-water 
ratio with water velocity parameter. The data follow smooth 
curves except, for the 6.29 cp. oil, water velocity 
0.1 Here again, the pressure gradients were higher 
than expected and can perhaps explained the same way 
for the 6.29 cp. oil. 


The smooth curves were used produce Figure which 
pressure gradient plotted against input water-oil ratio with 
oil velocity parameter. The flow pattern loci from Figure 
are superimposed. This figure shows the effect increasing 
the water content the flowing mixture constant oil flow 
rate. With this more viscous oil the pressure gradient reduction 
obtained the introduction water more pronounced. 
The maximum pressure gradient reduction factor noted was 
6.0 for oil velocity 0.5 ft./sec. input ratio 0.65 
(or water percentage 39.4% based the total liquid flow). 
Again high input water-oil ratios the observed pressure drops 
were reasonably close the calculated values assuming that 
the liquid flowing was all water. The pressure gradient was 
reduced most markedly the concentric and oil-slugs-in-water 
flow pattern regimes. smaller reduction was evident the 
water-drops-in-oil regime and possibly accounted for 
reduction the effective viscosity the water-in-oil emulsion 
high water contents below the viscosity the pure oil. 

Figures similar the figures for the less viscous oils are 
also presented for the 65.0 cp. the plot 
pressure gradient against input oil-water ratio with water veloc- 
parameter. Again the data fit smooth curves and with 
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Figure 18—Pressure gradients for the oil viscosity 65.0 
centipoise function the input water-oil ratio and oil 
velocity. Points are from the curves Figure 17. 


this, the most viscous oil, the high pressure gradients obtained 
when the oil formed the continuous phase are evident high 
input oil-water ratios. For given water velocity the pressure 
gradient increased with increasing input oil-water ratio with 
rather sharp increase the point inversion between the 
flow patterns with water the continuous phase and the flow 
patterns with oil the continuous phase. 


The practical significance the addition water the 
65.0 cp. oil better illustrated Figure where the pressure 
gradient plotted against the input water-oil ratio with oil 
velocity parameter. The flow pattern loci from Figure 
are superimposed. The curves all have the same general shape: 
gradual increase pressure gradient with increase water 
content through the flow regimes which oil the continuous 
phase, pronounced reduction pressure gradient when the 
inversion flow pattern with water the continuous phase 
takes place, and general increase pressure gradient with 
the continued increase input water-oil ratio. 

The maximum pressure gradient reduction factor for the 
65.0 cp. oil approximately and occurs for oil velocity 
0.3 ft./sec. input ratio 1.0. The superimposed flow 
pattern regimes indicate that the greatest pressure gradient 
reduction occurs the concentric oil-in-water and oil-slugs-in- 
water regimes. Again the observed pressure gradients high 
input water-oil ratios compare very closely with the pressure 
gradients calculated the assumption that water alone occupies 
the pipe. 

The maximum pressure gradient reduction factors obtained 
experimentally, i.e., 1.7 for the 6.29 cp. oil, 6.0 for the 16.8 
oil and 10.0 for the 65.0 cp. oil, may compared with the 
maximum reduction factors predicted the theoretical analysis 
Russell and for laminar concentric flow circular 
pipe. Russell and Charles predicted the maximum reduction 
factor 


(2u. = 


\O; - Mh. SO“. « in water 
/ | | — 


Substitution the viscosity values this expression gives 
maximum reduction factors approximately 3.8 for the 6.29 cp. 
oil, 9.6 for the 16.8 cp. oil and 36.4 for the 65.0 cp. oil. The 
experimental quantities are, therefore, not large the theoret- 
ical values. The difference explained noting that the flow 
pattern for which the experimental reduction factors were ob- 
tained was closer slug flow than concentric and, although 
the flow the water was probably laminar the annulus 
between the oil slugs and the pipe wall, eddies existed the 
water phase between the slugs. Nevertheless, the maximum 
pressure gradient reduction factors determined experimentally 
are considerable and show definite increase with oil viscosity. 


Conclusions 


The following flow patterns were noted the horizontal 
flow equal density oil and water mixtures: 


(a) water-drops-in-oil. 

(b) concentric oil-in-water. 
(c) oil-slugs-in-water. 

(d) oil-bubbles-in-water. 
oil-drops-in-water. 


The character the flow patterns and their relationship the 
oil and water flow rates were found largely independent 
the viscosity for oils viscosities 6.29, 16.8 and 65.0 cp. 
although the most viscous oil showed dissimilar behaviour 
high oil-water input ratios. This was probably due different 
oil-water interfacial properties. 


Holdup ratios were, general, found greater than 
unity when the water was the continuous phase contact with 
the pipe wall, less than unity when oil was in.contact with the 
pipe wall and close unity high flow velocities. 


The addition increasing amounts water oil 
which originally laminar flow (Reynolds number 
probably not more than 1500) lowers the pressure gradient 
minimum, after which the addition more water increases 
the pressure gradient and, with sufficient water, the pressure 
gradient exceeds the pressure gradient for the oil flowing alone. 


Nomenclature 


diameter pipeline, ft. 

Fanning friction factor, dimensionless 

dimensional conversion factor sec.? 

holdup ratio, the ratio the input oil-water ratio 
the situ oil-water ratio, dimensionless 

length pipeline, ft. 

input oil-water volume ratio, 

input water-oil volume ratio, 

Reynolds number, dimensionless. 

superficial liquid velocity, ft./sec. 

Subscripts 

force 

mass 

oil 

water. 
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theoretical study the countercurrent flow 
dispersed phase continuous phase resulted 
the following equation relating holdup the flow 
rates the two phases. 


For packed towers this equation has been used 
the following form: 


Values A’, B’, and for three liquid pairs, and two 
packing types are reported, using new holdup data 
determined this research. 

Graphical correlations are shown and values 
are reported for some the data Gayler and 
Pratt. 

For spray towers the equation suggested corre- 


lation 


has been tested using data previously published 
one the authors. 


0.2 


this paper defined the fraction the total 
liquid flowing packed spray type liquid-liquid 
extraction column occupied the dispersed phase. 


The importance holdup related the mass 


which the product the transfer rate coefficient 
and the interfacial area per unit volume. For given drop 
size the latter proporticnal the holdup. 
would impossible obtain values the transfer rate 
without knowledge the holdup and drop 
characteristics. 


Previous Work 


Holdup data the literature can classified into two cate- 
gories depending the manner which was obtained. 
the drainage method, measurements were made closing off 
flows the column, draining the column and measuring each 
the displacement method, all streams and from 
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the column are closed off, and then continuous phase added 
until the interface has been restored its initial 

Holdup has also been estimated observing the velocity 
individual droplets and calculating the number drops held 

Gayler and extensive study, observed that the 
holdup measured the displacement method, which they termed 
the normal holdup, was than the holdup measured drainage. 
The difference, which has been called permanent holdup, was 
presumably caused droplets trapped the packing. Gayler, 
Roberts, and have proposed the following equation 
relating holdup and superficial velocities below flooding 


Wicks and investigated the holdup toluene 
water using rather wide range column diameters and 
packing sizes. interesting note that they found that the 
reproducibility holdup results was improved greatly blow- 
ing the column with air after dumping the packing. This has 
been found the present research also. These authors proposed 
empirical equation the type 


Here the total holdup, the sum the normal and permanent 
holdup defined above. this equation the first term the 
right side predominates and accordingly there some difficulty 


Correlating Equation Proposed This Work 

(19) the paper Sakiadis and 
relating the holdup the dispersed phase the physical prop- 
erties the system and the flow rates both phases will 
used starting point for the derivation the holdup correla- 
tion. similar starting equation presented Dell and 
This equation was derived applying Bernoulli’s Theorem 
each phase separately and may written follows: 


Sakiadis found the value the order 0.2 for 
flooding conditions. However, conditions other than flooding 
difficult determine the value theoretically. Con- 
sequently was decided rearrange Equation (3) the 
following form: 


Mc" Cppp'-* Up?-* Mp" bp'** 
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Figure 1—Flowsheet liquid-liquid extraction apparatus. Heavy phase dispersed. 


Dividing through gives 


The groups containing the physical properties the phases 
were considered constant for any one system and type 
packing. Equation (5) may then written more simply 


( ( 1 ~ z) 


This equation has three constants for any system and packing: 
obtain these, the equation suggests that values 
assumed until correlating line obtained when 


this paper this equation will examined using spray 
tower data previously packed tower data obtained 
this and packed tower data Gayler and Pratt. 


Apparatus and its Operation 


flow diagram the apparatus used for holdup measure- 
ments packed towers shown Figure The arrangement 
essentially the same that used the spray tower studies 
whose results are also used this paper. The major difference 
the replacement the swinging loop needle valve and 
flowrator for more precise the dispersed phase flow 
leaving the column. 

Figure shows details the tower. The column proper 
consisted 2-ft. length Pyrex pipe i.d. The spray 
nozzle was made cluster eight tubes arranged 


circle diameter. The packing support was con- 
structed wide strips type 302 stainless steel sheet 
mounted two rods and separated porcelain 
Raschig rings spacers. 


The packings used were unglazed porcelain Raschig 
rings and Intalox saddles supplied the United States Stone- 
ware Company, Akron, Ohio. 

For the packed tower runs, holdup was measured using 
combination the displacement and drainage methods 
follows: 


The two phases were brought the operating temperature 
25°C. the use the heat exchangers shown Figure 
and each liquid was saturated with the other before run was 
started. The column was filled with the continuous phase and 
after few minutes circulation the dispersed phase was 
introduced the top the column. The column was allowed 
reach steady state and operate this condition for 
min., the dispersed phase being removed the bottom 
rate maintaining constant interface level the settling chamber. 

All flows and from the column were shut off simultane- 
ously. The dispersed phase the packing was allowed 
settle and was then drained and measured graduated cylinder 
while equivalent valume continuous phase was added 
the column. This measured the normal freely falling holdup. 


The column was then drained and washed several times with 
continuous phase remove any the permanent holdup trapped 
the packing interstices. 

The sum the free falling and permanent holdups yield the 
total holdup. 


Presentation Holdup Data 


(a) Packed tower data 


The packed tower holdup data reported here were measured 
during the course studies transient behavior liquid-liquid 
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Figure 2—Extraction column. Heavy phase dispersed. 
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Figure holdup data. System: water; 
Packing: Tutalox saddles 0.756. 
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Figure 5—System: CCl, water; Packing: Raschig 
rings. 
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Figure 6—Final holdup correlation. System: CCI,—water; 
Packing: Raschig rings 0.587. 
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rings. (Data Gayler and Pratt). 
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Figure Butyl acetate-water. Packing: Raschig 
rings. (Data Gayler and Pratt). 


extraction packed towers. Accordingly the intention was 
obtain precise data over rather small range flows. 


Previous authors found that the method packing column 
had appreciable effect the reproduceability holdup data; 
general, data deviated much 20% after repacking. 
Wicks and settled the packing passing air 
through randomly packed column full water rate 
sufficient cause the packing surge. this study this 
method repacking resulted columns that were reproduce- 
able 13%. 
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Figure carbitol-water. Packing: Raschig 
rings. (Data Gayler and Pratt). 


Typical holdup data for carbon tetrachloride dispersed 
water are shown Figure The free falling and total holdup 
both appear linear functions the disperse phase rate 
given continuous phase rate the loading point. This 
agreement with the findings previous 
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Figure 10—4-in. i.d. spray tower data. System: Ethylene dichloride-water. 
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permanent holdup data, which the difference between 
the holdups shown this graph, increases with increasing 
dispersed phase rate well with increasing continuous 

hase rate. 

The holdup the Intalox saddle packing was found 
less than that the Raschig ring packing. This agreement 
with other who found that holdup decreased 
with increase the fraction void packing size. 

total holdup data below flooding for three systems and 
two packing shapes were correlated suggested Equation 
(6). indicated Figure when Equation (6) was 
taken 0.5, straight lines differing slope but with common 
intercept could drawn. Assuming that the slope would vary 
exponentially the continuous phase rate, Figure was drawn. 
This led the following modified form Equation (6). 


Up! 5 Up'5 (; z) ( ) 


Based this equation the total holdup correlation for carbon 
tetrachloride water shown Figure 

Table shows values A’, and Equation (7). 
would appear the may function the type packing. 
would expected that and are functions the physical 
properties the systems but there insufficient data this 
study find the relationship, any. 


(b) Packed tower data Gayler and Pratt (2) 

For comparison purposes and test the generality the 
correlating equations proposed this paper the data for 
Raschig rings Gayler and were examined. Plots 
such those shown Figures and this paper were 
first prepared, and the final correlating graphs shown Figures 
and were produced. 

seen that compared with Figure considerable curvature 
was obtained, perhaps due use data covering higher range 
holdup values and also that less satisfactory correlation 
the data for carbitol-water than for the other two systems 
was obtained. 


(c) Spray tower data 

Holdup data for the system ethylene dichloride-water, deter- 
mined i.d. spray column the Elgin have 
been correlated the method proposed this investigation. 
this case was taken have the value 0.2 proposed 

indicated Figure 10, four different curves were ob- 
tained, rather than straight lines. The reason for the curvature 
may lie the wider range flow conditions covered, these 
data extending the rejection definition The 
four curves represent four different nozzles, producing different 
drop sizes and distributions. 

the value the ordinate appears the abscissa and 
one might have expected greater tendency for the data 
yield graph slope unity. The fact that such slope was not 
obtained indicates that the plotting technique not false one 
and that the and values Equation (7) are significant. 
However, this fact, along with the high powers and 
that appear the correlations, will limit the usefulness the 
correlating technique. 


Conclusion 


this paper theoretical approach holdup below flooding 
has been proposed, resulting Equation (6). practice this 
equation has been modified Equation (7) which must 
regarded empirical. 

Some new data for packed towers has been presented the 
proposed correlation. Previously reported spray tower data 
have been handled the method also. 
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TABLE 
HOLDUP CORRELATION 


x8 A Ug x B 


1/2 in. Raschig ring packing 


System 
Carbon tetrachloride—water 0.023 0.000031 0.61 
Ethylene dichloride —water 0.0093 0.000152 0.64 
Carbon tetrachloride 
naphtha—water 0.0114 0.65 


1/2 in. Intalox saddle packing 


System 
Carbon tetrachloride—water 0.343 0.0000026 0.12 
Ethylene dichloride —water 0.248 0.12 
Carbon tetrachloride 
naphtha—water 0.244 0.12 


Nomenclature 
Roman 
Fractional voidage packing 
Gravitational constant ft./hr.? 
Transfer coefficient ft./hr. 
Periphery terms for the continuous 

dispersed phases, defined Sakiadis 

Superficial velocity disperse phase ft./hr. 

Superficial velocity continuous phase ft./hr. 
Acharacteristic velocity ft./hr. 
Fractional holdup 
Greek 
Viscosity disperse phase ft. 
Viscosity continuous phase ft. 
Coefficients 

Exponents 
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Approximate Solutions Conduction 
Heat through Non-Homogeneous 


CHI TIEN? 


comparing the respective difference equations 
heat conduction through homogeneous and non- 
homogeneous media, approximate analogy 
established. This analogy enables the finding the 
solution conduction through non-homogeneous 
medium terms the known solution for homo- 
geneous cases. illustrative example given 
demonstrate the working details. 


study the mathematical solutions heat conduction 
has been subject interest for many years. Numerous 
solutions corresponding various boundary conditions are 
available and these have been summarized the excellent 
work Carslaw and However, most these solu- 
tions apply only cases where the conduction takes place 
through medium. the case non-homogeneous 
medium where the physical properties are functions location 
although there difficulty obtaining the solution the 
strict mathematical sense, the solutions are often expressed 
terms infinite integral. The evaluation these integrals 
often requires elaborate and tedious computation. This has 
been proved impractical some the engineering applications. 

The object this investigation propose simple method 
which can used obtain approximate solutions for such 
problem. This method based the analogy between con- 
duction homogeneous medium and that non-homo- 
geneous one, which established considering the respective 
difference equations. 


Mathematical Development 


For the case one-dimensional conduction non-homo- 
medium where the physical properties are functions 
location, the conduction equation Cartesian coordinates can 
written as: 


The corresponding difference equation can given as: 


Tings — Tin 
TS 


Axm m 


Where and denote the magnitude the position and 
time increments respectively; the subscript indicates the value 
the mth position increment while the superscript designates 
the value the jth time increment. should noted that 
although the time increment being kept constant, the magnitude 
iManuscript received September accepted October 1960. 
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the succeeding position increment not the same but varies 
with given relationship given later. Equation (2) can 
written as: 

Axm Axm+i( Axm +Axmsi) 

2 Axa + Rm—1/2 Axm+i) Tne] 

— 


[Rm+1/2 Axn Tingitkm—12 Axm+1 


(3) 


For sufficiently small increments, the following approxima- 
tion believed applicable: 


Axe = Rm—1/2 Axm+1 = 1/2Rm( Axm + Axn+1) (4) 
Substituting Equations (4) into (3), results: 


Furthermore, Equation (5) would simplified the following 
form: 


the following condition satisfied 


Axn+1 Axn 


where 


For the case homogeneous medium, the corresponding differ- 
ence equation known as: 


and 
=i 


Equations (6) and (10) are identical form except that 
the latter case (homogeneous), the magnitude constant 
while the case varies depending the 
particular location and the recurrent formula given Equation 
(7). This identity functional forms leads the conclusion 
that with proper choice position increments, the temperature 
history one point the homogeneous medium and that 
another the medium are the same provided 
there exists similar initial and boundary condition these 
two cases. initial condition, does not necessarily 
mean that the two initial conditions are identical; implies 
rather that the temperature the corresponding points are the 
same. This conclusion can also stated follows: With 


The Canadian Journal Chemical Engineering, February, 1961 


similar 
arbitrary 


increment 


where 
Ax; = 
With this 
heat 
the 
available. 
the specia 
similar 
tion 
obtain 
geneous 
case shot 
available 
the 
cross 
choice 
since the 
totally 
regard, 
least 
limiting 
taken 

Vario 
ful. How 
linear 
relationsh 


Assun 
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= = 
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2 
The Can 


= 


similar initial and boundary conditions, the temperature history 
K.(Ax), for homogeneous medium, where any 
arbitrary integer and (Ax), designates the convenient space 


increment, and that for non-homogeneous medium, 
= 

where Ax; and are related Equation (7) and Ar, 
(Ax,) with the thermal diffusivity the 
homogeneous medium and the thermal diffusivity the 
surface the non-homogeneous medium, will the same. 
With this conclusion, becomes possible estimate the solution 
heat conduction through non-homogeneous medium provided 
the similar solution for the case homogeneous medium 
available. Furthermore, although this conclusion based 
the special case unidirectional and with Cartesian coordinates, 
similar results may obtained for other cases too. The applica- 
tion this method then becomes fairly simple. order 
obtain the temperature distribution throughout 
medium with its physical properties given function 
position, the corresponding solution for the homogeneous 
case should obtained first. This information 
available the literature. With the information, one can obtain 
the temperature history several selected points the non- 
homogeneous medium with the use the method stated before. 
cross-plotting this temperature history, the temperature 
distribution any time can obtained. obvious that the 
choice the initial position increment, Ax,, can very critical 
since the validity the approximation given Equation (4) 
totally dependent the magnitude the increment. this 
regard, suggested that results should obtained using 
least two different values for the initial increment and the 
limiting value obtained through extrapolation should 
taken the actual value. 

Various attempts have been made express this relationship 
explicit form. Unfortunately this has not been success- 
for the special case where thermal diffusivity 
linear function the depth the medium, more direct 
relationship can obtained. This will shown follows: 

Assume that the thermal diffusivity this particular medium 
has the following form: 


The corresponding locations the homogeneous and non-homo- 
geneous media, where the temperature histories are similar, are: 


Substituting Equation (11) into (16) for have: 


Ax. = 


G2 At a[l + a(Ax Ax:)| 


Ax; 2 At = 
Ax 
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general, one obtains: 


Substituting Equations (19a) and (19b) into (14), obtain: 
N-2 


s=1 


N(2a, At) for even 


N-1 


2 


now seek the limit the ratio, for given value 


N 
Ax; 


li P t=1 

N — © 

constant for even number 


N — 


4a, 


or 


tl 
4+ 


if = Gy... 


Solving the quadratic equation Equation (21a), have: 
a 


With this relationship established, one can obtain the temper- 
ature distribution non-homogeneous medium where 
the solution for the case where a,. illustrative example 
will given the following section. 


EXACT 


APPROXIMATE 
SOLUTION 


= 


EXACT SOLUTION 
APPROXIMATE SOLUTION 


2.0 3.0 


Figure 1—Comparisons between the exact and approximate solutions for two sets conditions. 


Example 


For illustration, consider the following problem: semi- 
infinite solid where initially zero temperature 
imposed the surface. The exact solution for the temperature 
distribution given Carslaw and is: 


Numerical values are available Reference (2) 


the approximate method applied, first for similar 
problem homogeneous medium where the tempera- 
ture distribution given by: 


(24) 


Substituting Equation (22) for the Equation (24), one 
obtains the approximate solution as: 


T, av at 


Comparison between Equation (23) and (25) shown 
Figure For both cases the approximate solution gives 
higher numerical value. conceded that the degree accuracy 
not good expected (about 10% for the range 
.01/—0.9), the time saved computation quite significant. 

For some preliminary engineering calculations, the use 
this method may more advantageous. 


Nomenclature 
Constant 
Heat capacity 


Superscript designating the value the jth time incre- 
ment. 


Thermal conductivity 

Subscript indicating the value the mth position incre- 
ment. 

Positive integer 

Tempetature 

Time 

Distance 

Thermal diffusivity 

Time increment 

Position increment 

Density 
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NOTE THE EDITOR 


The Extension Couette Flow Solution Fluid 


CHI TIEN 


Essex College, Assumption University Windsor, 
Windsor, Ont. 


purpose this note show that the solution 
Couette flow problem laminar heat transfer for 
Newtonian fluid can also applied the case non-Newtonian 
system, provided that the relevant parameters are properly 
defined. 

Couette flow, refer the case flow between two 
parallel plates, which one stationary, while the other 
moving constant velocity prescribed thermal condition 
imposed these two plates and the question obtain the 
temperature distribution. For the treatment this problem 
with Newtonian fluid, one can refer standard heat transfer 
textbook, i.e., 

The two-dimensional energy equation incompressible, 
constant property fluid given as: 


For Couette flow, when pressure drop along x-direction 
between these two plates. 

The last term represent the heat dissipation due friction. 
For case like this, given by: 


Although there are general postulates relating the stress 
and strain tensors for non-Newtonian fluid, for large class 
material, power model relationship can used, i.e., 


Substituting (2) and (3) into (1) together with the velocity 
distribution, obtain: 


de l+n l+n 


the two plates are maintained constant temperatures 
such as: 
T=T, 
T=T aty=h 


Integrating Equation (4) twice, the temperature distribution 
given as: 


Rewrite the dimensionless coefficient the second term as: 


The second group the right known Eckert number. 
The first group can written as: 
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Since 
pin 
for fluid. 


one assumes that the same reationship exists between the 
various dimensionless groups for non-Newtonian fluid that 
Newtonian fluid, such Np, Equation (5) can 
put the following form: 


This expression identical that the Newtonian fluid. 

Another case worthwhile considered the one 
which one plate insulated while the other maintained 
constant temperature such as: 


The temperature corresponding this condition can 


shown as: 


The temperature difference between the bottom and top 
plate, result adiabatic heating given as: 


=T,-—TI,= = = Np,—.. 
the recovery factor defined as: 

= 


These two expressions are again identical form with those 
Newtonian fluid. 


Nomenclature 
Width between plated ft. 
Parameter the Power Law model for non-Newtonian 
fluid 
Exponent the Power Law model for non-Newtonian 
fluid, dimensionless 
Echert number defined U?/C,AT 
Np, Peclet number defined C,Uph/k 
Np, Prandtl number defined Np, 
Reynold number defined 
Temperature °F. 
Velocity along x-direction ft./hr. 
Velocity the lower plate ft./hr. 
Velocity along y-direction ft./hr. 
Velocity vector 
Density 
Stress tensor 
Reference 


(1) Schlichting, H., “Boundary Layer Theory”, Pergamon 
Press (1955). 


INDUSTRIAL 


The Explosive Limits and Flammability 
Xanthate Dusts Air’ 


The minimum concentrations xanthate dusts 
air which will form explosive mixtures have been 
determined. Comparison with other dusts considered 
hazardous industrially indicates that xanthates 
the lower alcohols should classified readily 
explosive. The explosive limit increases rapidly with 
particle size; dusts coarser than 150 mesh are rela- 
tively safe. Xanthate dusts are easily ignited and have 
very low ignition temperatures. Under conditions 
such that fine wood flour ignites 540°C. the 
xanthates will ignite between 230°C. and 290°C. 


acids are the O-esters dithiocarbonic acids. The 
sodium and potassium salts xanthic acids derived from 
the lower aliphatic alcohols are water soluble and stable and 
find extensive use froth flotation. They are usually manufac- 
tured reacting sodium potassium hydroxide with exess 
the appropriate alcohol and carbon disulphide. Purity the 
commercial product may range from with dithiocarbon- 
ates and thiosulphates the main impurities. The free xanthic 
acids are unstable but their alkali metal salts may handled 
relatively safely either pellet powder form. generally 
that they should kept cool and dry prevent 
excessive decomposition, and cases spontaneous combustion 


have been reported. 


reference the flammability explosive tendencies 
xanthates reported the literature, but since their decomposi- 
tion both heat and moisture known regenerate much 
the carbon disulphide and alcohol might expected that they 
will burn and their dusts explode fairly readily. 


The hazard fire from flammable dusts usually less sig- 
nificant than the explosion hazard and occurs only when layers 
are allowed accumulate. The temperature required ignite 
such establishes the degree hazard and some cases 
low enough that low temperature sources ignition such 
overheated bearings and motors, light bulbs and steam lines 
must considered 

The hazard dust explosions given manufacturing 
operation depends many factors. The minimum explosive 
concentration dust air depends primarily the material 
constituting the cloud and the particle size distribution it. 
addition, the type ignition source and its temperature will 
affect this minimum concentration, will such factors the 
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relative humidity and moisture content and the presence 
absence air currents. Generally speaking, industrial interest 
centered that minimum concentration which will explode 
under optimum conditions all other factors. For this reason 
most practical test dusts under conditions low moisture 
content and relative humidity still air using high energy 
source ignition. Finer and finer fractions the dust should 
tested until little further effect from decreasing particle 
size observed. 

Under certain conditions the assessment hazard and design 
prevention techniques may require knowledge character- 
istics the dust which are purely academic interest other 
cases. The minimum spark energy necessary ignite dust 
clouds under various conditions and the pressures developed 
the explosion are known for many Where inert 
atmospheres are maintained process equipment the 
minimum oxygen content which dust explosion can 
induced will important. 


Certain broad limits can set the conditions under 
which dust cloud likely hazardous. Most flammable 
substances which are not pyrophoric true explosives are not 
hazardous the particles are all nearly all coarser than will 
pass 200 mesh screen. Large particles settle quickly still 
air and are difficult ignite suspension. The lowest concen- 
tration fine dust which will explode usually the range 
100 mg. per liter air, but some resins will explode concen- 
trations low mg. per liter. Very few substances will 
explode atmospheres containing less than 10% oxygen. 


establishing the potential from dust explosion and 
ignition settled dust for specific compound the method 
testing must carefully selected. Ideally, infinitely large 
cloud uniform composition still air required establish 
lower explosive limit for dusts. Practically, small chamber 
about one liter used and the dust dispersed sudden puff 
air. comparison with other substances known explosive 

tendencies value can obtained such equipment that will 
close the true lower explosive limit, but care must 
taken ensure good dispersion the dust and adequate energy 
and temperature the igniting source. The detection explo- 
sions near the explosive limit sometimes doubtful and some- 
what subjective. 


Ignition temperatures layers dust may subject 
similar experimental variations. The rate heating and the 
length time that the sample held the test temperature 
are important. Air currents and the materials contact with 
the dust will also affect the temperature which ignition occurs. 
Again, comparison with materials generally known hazard 
tested under the same conditions seems advisable. 
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Figure 1—Particle size distribution commercial samples 
xanthates. 


Experimental and Results 


The explosion chamber used for these tests was described 
earlier weighed sample the dust tested 
was dispersed metered puff compressed air one 
liter transparent plastic chamber. 5000 volt a.c. spark was 
used induce explosions. positive test was recorded the 
paper covering the top the chamber was ignited burst. 
Ignition temperatures were determined placing mg. the 
sample stainless steel plate 8in.) the immediate 
vicinity thermocouple silver soldered the plate. burner 
was placed under the other end the plate and regulated 
that the temperature rise the thermocouple was 80°C. per 
minute. The output from the thermocouple was recorded. The 
ignition was readily observed the chart sharp rise 
temperature. 


Samples dust were equilibrated various relative humid- 
ities determine moisture content and its effect explosive 
limits. Relative humidities were kept constant saturated 
solutions suitable salts the humidity chambers. 


Separation the samples into fractions known particle 
size was carried out Roller Particle size Analyzer* modi- 
fied for higher air flows. 


Samples four xanthates commercial grade were separ- 
ated into six fractions. The particle size distribution all four 
samples was similar (Figure 1). noted that from 
70% each the samples tested was the size range which 
could expected cause dust explosion. The samples 
received were essentially anhydrous and lost less their 
weight prolonged vacuum drying room temperature. 

preliminary test, the minimum explosive concentration 
air the fraction passing 250 mesh screen was determined 
for three these xanthates. Results were somewhat variable 
because the dusts failed disperse well the test chamber 
but values the range from mg. per liter for sodium 
xanthate mg. per liter for sodium isopropyl xanthate were 
recorded. tests were conducted time when the 
ambient relative humidity was very low (less than 10%) and 
later could only duplicated samples dried under vacuum. 


Potassium xanthate agglomerated such extent 
that could not screened successfully, and the other xanthates 
were nearly bad, further fractions were obtained air 
separation. 

Figure shows the effect particle size the explosive 
limit vacuum dried fractions. Even the fraction finer than 
microns did not have explosive limit lower than the nominally 
much coarser 250 mesh screened tested, 
possibly indicating some degradation these compounds the 
course air separation. expected, the explosive limit rises 


*American Instrument Company, Silver Springs, Md., U.S.A. 
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Figure 2—The effect particle size the minimum 
explosive limit xanthate dusts air. 


sharply for coarser particles, and beyond the 100 micron (150 
mesh) level the hazard represented small. 


The xanthates take considerable moisture but quite slowly. 
Potassium xanthate does not absorb appreciable moisture 
relative humidities less than 50% but the other three absorb 
from relative humidities 12% 3). Sodium 
xanthate comes the same moisture content all levels 
humidity tested whether the sample originally had more less 
than the equilibrium amount, but sodium ethyl xanthate and 
sodium xanthate exhibit hysteresis the lowest levels 
relative humidity tested (Figure 3). possible that this 
results from the formation hydrate the higher humidities 
which, once formed, stable the lower humidity. 
jection the descending moisture content relative humidity 
relation shown for the three sodium xanthates 
gives intercept case very close the theoretical for 
the moisture content the monohy drate. 

Moisture absorption the highest relative humidity tested 
(77%) continues very slowly after the hours used determine 
moisture contents these tests, and after some 
weeks sodium ethy] xanthate will dissolve the water absorbed. 
This slow absorption accompanied decomposition the 
compound, probably through drolysis carbonate and tri- 
thiocarbonate. Potassium xanthate absorbs moisture very 
slowly and after one month 77% relative humidity contains 
less than 15% moisture. 

Because the moisture absorbed, the explosive limits 
these three xanthates are very much higher equilibration 
even low relative humidities (from higher moisture contents) 
than the limits given for the anhydrous compounds (Table 1). 
Since the effect moisture typical and maximum hazards 
are most concern single example shown. Equilibration 
slow (Figure and even thin layers not complete for 
many hours room temperature that the sample may exhibit 
explosive properties characteristic the anhydrous compound 
for some time after exposure ambient conditions, particularly 
relative humidities the range where hysteresis may occur. 


TABLE 


MOISTURE THE MINIMUM EXPLOSIVE 
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Figure content xanthates after exposure 
atmospheres various relative humidities for hours. 


TABLE 


IGNITION TEMPERATURES XANTHATES 
AND OTHER SUBSTANCES 


Compound Ignition Temperature 

Ethyl Xanthate 245 
Isopropyl Xanthate 270 
Xanthate 260 
Amyl Xanthate 280 
Sulphur 250 
Wood Flour 550 
Starch >800 


The ignition temperatures the four xanthates tested are 
given Table along with ignition temperatures for fine wood 
flour and sulphur, determined the same way. Starch was also 
tested but failed ignite Since wood flour and 
sulphur have ignition temperatures among the lowest recorded 
for industrial dusts, the xanthates can said have very low 
ignition temperatures. The ignition temperature found this 
method for wood flour higher than reported the 
which may result content the variety wood 
used, but the ignition temperature for sulphur agrees with 
reported values. Moisture content did not affect the ignition 
temperature the xanthates any large extent. 


Summary 

Commercial samples four xanthates contained 70% 
material fine enough constitute potential dust explosion 
hazard. Materials received were anhydrous but three appeared 
form rather unstable monohydrates when exposed atmo- 


77% REL. HUMIDITY 

2 
XANTHATE 
ETHYL 
sec-BUTYL 

° 12 24 

TIME, hours 


Figure 4—The rate uptake moisture two xanthates 
exposed atmosphere with relative humidity 77%, 


spheres with relative humidity above 50%. Two the 
xanthates exhibited hysteresis the lower regions relative 
humidity but sodium xanthate had the same equilibrium 
moisture content relative humidity whether equilibrium 
was approached from higher lower moisture levels. 

Explosive limits for the finest fractions the anhydrous 
xanthates dispersed air were the range mg. per 
liter. The samples containing moisture equivalent the mono- 
hydrate had explosive limits about 100% higher and were more 
difficult disperse the test chamber. 

The ignition temperatures for these four xanthates were 
the range 280°C., comparable with such hazardous 
materials sulphur. The method used determine ignition 
temperatures would not ignite starch 800°C. and would 
barely ignite fine wood flour (550°C.) 

the anhydrous state these xanthates represent distinct 
hazard with respect dust explosions, particularly since they 
contain large amounts fine material and absorb appreciable 
amounts moisture only very slowly. Their ignition tempera- 
tures are low enough rank xanthates one the more 
flammable commercial chemicals. 
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Method Measuring the Surface 
Area Granular Material’ 


method described which permits measure- 
ment the exterior surface areas objects down 
about 0.1 sq. in. area. this method, the object 
coated with adhesive and single layer small 
metal particles. The change weight gives measure 
the external the object. 


Systematic variations the density the metal 
particle coating are discussed, and appropriate cor- 
rection factors are given. 


This method intended for use research that 
involves objects whose exterior surface areas have 
effect the process being investigated. Examples 
are the drying agricultural grains and foodstuffs, 
the comminution ores and chemicals, and the 
hydraulic and pneumatic transport materials. 


The precision the method adequate, the root- 
mean-square error being about 1%. 


measurement the surface areas irregular objects 
has occupied the attentions number investigators. 
Many the methods have been reviewed Dalla 
and Orr and Dalla 


The term “surface area” requires careful definition and its 
numerical value may depend upon the method measurement. 
For example, necessary distinguish between exterior 
surface area and one associated with pores and fissures within 
the body; the technique measurement would depend upon 
which area” was required. 


The method that the subject this paper consists 
attaching single layer metal particles 
the surface the object which the area determined, 
(Fi igure and measuring the attendant change weight. There 
obvious upper limit the size object, and the method 
applicable for surface areas down about 0.1 sq. in. Such 
objects may tested singly groups. 


The metal particles adhere the exterior, visible surface 
the object, and the area this surface that measured. 
This type surface, the above mentioned range surface 
areas, not readily measured other techniques, such 
adsorption permeability. The adsorption technique applic- 
able much smaller particles and, larger total areas, and includes 
the areas pores and fissures. The permeability method 
indirect and depends upon the packing and porosity bulk 


material. 


\Manuscript received January 6; accepted October 21, 1960. 
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When formulating expression describe the behavior 
particles engineering process, frequently necessary 
include characterizing phy sical dimensions the formula. The 
technique discussed here may useful investigaticns 
processes which the exterior surface area plays part. This 
area may enter the formula directly, may used derive 
shape factor characterize the particles. Examples are: the 
motion particles gases and liquids, the rates drying 
granular materials, the rates chemical reaction particle 
surfaces, and the crushing and grinding ores and chemical 
substances. 


Method Coating 


the selection the powder, several factors were con- 
sidered. The particle size must small compared with the 
object measured, and indeed should small compared 
the smallest radii curvature that commonly occur the 
object. The particles should have high density provide 
easily measured change weight and, finally, they should not 
stick together. The powder actually selected met these require- 
ments satisfactorily. was pure nickel powder, with particles 
fairly regular form. quantity was procured* and sieved 
obtain ‘particles the sieve size range —270 +325 (US. 
Standard Sieve) (2.1 1.7 mils diameter). 


The critical problem this method that attaching the 
particles single layer uniform and reproducible density. 
Success depends primarily the selection and proper application 
suitable adhesive. Experience showed that the thickness 
the adhesive layer had considerably smaller than the 
diameter the powder particle order prevent multiple 
layering. The best results were obtained using slow-drying 

uniform layer adhesive about 0.0002 in. thick was 
placed the object first dipping the varnish, then 
shaking vigorously about 500 gm. mesh sand 
until the sand did not adhere it. (10 seconds usually 
sufficed). Since dry sand removed too much the varnish, 
was necessary prepare advance thoroughly mixing 
gm. varnish with 500 gm. sand, and spreading the coated 
sand out dry air for about min. room temperature. 

The adhesive-coated objects were removed from the sand 
and were allowed dry air fine screen for 2-3 min. 
They were then placed vessel containing nickel powder 
and were shaken with for predetermined length time. 
The vessel was emptied onto wire screen over beaker 
recover the objects which were then rolled the screen 
tapped lightly remove unattached powder. The coatings 


pure nickel powder was kindly supplied Mr. Galley Sherritt- 
Gordon Mines Ltd., Toronto. 
*Berry Bros. Lionoil floor varnish was used in most of the tests. 


— 


INDUSTRIAL 


Figure 1—Nickel powder coating flax seed. 


proved quite securely attached and could handled with 
tweezers without damage. The objects were weighed 
analytical balance before coating with adhesive and after coating 
with nickel powder. 


Factors Affecting the Deposition Metal Powder 


During the development the method became apparent 
that number variables systematically influenced the weight 
powder attached per unit area. investigation was carried 
out isolate and study the most important effects arising from 
the two principal sources: (1) the properties the object 
being coated and (2) the experimental technique 
this purpose, variety regular objects was made materials 
such steel, glass, aluminum, Bakelite and Plexiglas, 
variety shapes and sizes. each case, the area and density 
the object was determined measurement. While many 
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Figure 2—The envelope density powder vs. reciprocal 
wire diameter. 


other objects were tried, the ones most used were among the 


(a) Spheres, diameters from 0.1 in. 0.8 in. 


(b) Cylindrical rods and wires, diameters from 0.003 in, 
0.25 in. and lengths from 0.25 in. in. 


(c) Square and rectangular bars with cross-sectional dimen- 
sions from 0.016 in. 0.2 in. and lengths from 0.9 in, 
in. 


(d) Flat plates, thicknesses from 0.02 in. 0.1 in. and 
other dimensions from 0.25 in. in. 


The various effects are discussed below. 


Time Shaking the Nickel Powder 


the coating process the surface density the powder 
increased rapidly the initial stages, and the rate increase 
dropped sharply available spaces were occupied the par- 


Envelope Effect 


For any surface with radius curvature much larger 
than the powder particle diameter this method estimates the 
actual surface area. may assumed that the particles occupy 
envelope area that stands distance equal one-half 
the particle diameter, from the object surface. such the 
case, the areas curved surfaces will overestimated the 
ratio (D+d)/D for rods and wires for spheres. 
This assumption may tested deducing the envelope thick- 
ness from experimental results and comparing with the 
particle radius 


let 


weight gain/unit actual surface area 


weight gain/unit envelope area 


shaker 


rotator 


specific volume object 


Figure 3—The envelope density powder vs. object 
specific volume for two methods coating. 
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series experiments was performed using rods and wires 
ranging diameter from 0.071 in. 0.0032 in. For these 
sizes the ratio (D+d)/D ranges from about 1.03 1.60, based 
powder particle diameter 0.0019 in. 

shows the graph vs. 1/D. The straight line 
was fitted the method least squares, giving the results: 

0.0095 in. 
121 


second series experiments performed with different 
wires gave: 
0.0011 in. 
123.8 


The agreement between the experimental values and 
the value d/2 sufficiently close support the use this 
concept the basis for practical correction factor. 


Density Effect 

The density the object being coated had small but 
significant effect upon the surface density the coating. 
investigate the influence object density the absence any 
shape size variation, identical steel cylinders were drilled 
axially with holes various sizes and filled with lead plastic. 
The result was series cylinders which, other than the density 
variation, had only slight variation the surface condition 
the ends. 

These objects, togenther with geometrically identical objects 
Plexiglas, Bakelite and aluminum, showed clearly that surface 
density coating increases with object density. 


Effect Method Coating 

The magnitude the density effect was influenced 
fourth variable, the method coating. methods 
applying the nickel powder were tried. each the object was 
placed vessel with amount the powder and was 
shaken. the first method, 100 ml. test-tube with 
nickel powder was shaken vigorously hand. the second, 
2000 ml. beaker containing 200-300 gm. powder was 
shaken bottle shaker that carried the beaker 180 r.p.m. 
horizontal circle in. diameter. the third, 1000 ml. 
containing about 300 gm. powder was rotated end over 
end r.p.m. 

Figure shows the results three runs each the 
variable-density objects the second and the third coatin 
methods. Both methods show the density effect, the shaker 
more strongly. 

When the shaker was used, the buoyancy effect the nickel 
powder kept the low density objects partly above the powder, 
and therefore less exposed packing influences. 

The variation density effect with coating method 
present precludes the application simple correction factor 
and can best accounted for including standard object 
having about the same density that the unknown object. 

Figure shows that, apart from the influence density 
effect, the results the two methods are very similar. 


Precision the Method 


find the precision the method, series six runs 
was made objects four different shapes (sphere, rod, 
bar, plate) and three different densities (steel, glass aluminum, 
plastic). (Figure 4). The third (rotator) coating method was 
used, and the results mg./sq. in. are given Table 
analysis variance shows that, while the density effect 
highly significant, the variations between shapes and between 
runs are not significant (their mean squares being quite close 
the error mean square). The analysis gave root-mean-square 
error 1.0%, 1.2 in. 
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Figure 4—Objects used obtain the results illustrated 
Figure (the rods and bars are about long). 


TABLE 


RANGES POWDER COATING DENSITY 
FROM EXPERIMENT 
mg./in.? 


Object material and density 


Object 

7.76 2.48 2.70 1.16 

Steel (Glass, Aluminum) Plexiglas 
Plate 120—121.5 —118.5 
Bar 121—123.5 119 —120 115.5—119 
Rod 122—124.5 117.5—118.5 
Sphere 123 117.5—119.5 117.5—119.5 
Summary 


method has been devised for estimating the external 
surface area for objects particular size range. This method 
was developed for research cereal grains but has applications 
research any particles similar size. For example, 
might used measure the surface area granular materials 
entering chemical process. 

this method the critical problem that placing 
suitable layer adhesive the object. The technique described 
here was satisfactory for the objects used this experiment. 
For other shapes, particularly those with negative curvatures, 
not entirely satisfactory and some other technique might 
prove better. 

Certain simple precautions may taken ensure the 
validity the measurements. Examination the coatings 
with stereoscopic microscope provides control the coating 
technique. The systematic variations due the nature the 
object are small, and can handled comparison method. 
That is, standard objects known area, roughly similar 
size, shape, and density the unknown objects, may included 
with them the coating and weighing techniques. experience 
gained, the amount control applied may reduced. 
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